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The  basic  performance  of  an  adaptive  antenna  system  is  influenced  some¬ 
what  independently  by  two  sub-systems:  the  antenna  and  the  adaptive  processor. 
Choice  of  the  antenna  type  (multiple  beam  antenna  or  phased  array)  and  the 
design  of  the  adaptive  processor  (which  is  used  to  control  the  weighting  at 
each  antenna  feed  port)  depends  strongly  on  the  specific  requirement  of  a 
particular  system.  To  date,  much  has  been  published  on  the  ideal  performance 
characteristics  of  adaptive  nulling  antenna  systems.  However,  little  has 
been  published  on  the  effects  of  hardware  component  imperfections  on  system 
performance.  To  characterize  these  effects  is  the  purpose  of  this  note.  We 
present  methods  for  categorizing,  analyzing  and  measuring  the  effects  of  a 
few  of  the  key  components  which  are  common  to  the  majority  of  adaptive  systems 
implemented  to  date.  Some  examples  we  will  discuss  are  those  effects  due  to 
antenna  type,  antenna  tolerance  errors,  channel  tracking  errors,  weighting 
circuits,  hard-limiters,  correlators  and  base-band  components.  We  analyze 
these  from  the  viewpoint  of  the  degree  of  cancellation  achievable  by  the 
adaptive  processor  and  also  their  impact  on  loop  dynamics  (where  appropriate) . 
In  order  to  illustrate  the  various  effects,  measured  performance  results 
obtained  using  an  experimental  Applebaum-Howell  type  adaptive  processor  will 
be  presented. 
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I.  INTRODUCTION 


Adaptive  antennas  are  capable  of  autonomously  changing  their  radiation 
pattern  characteristics  in  response  to  the  signals  incident  on  the  aperture. 
This  is  generally  accomplished  by  applying  amplitude  and  phase  control  (weight¬ 
ing)  at  the  output  of  each  of  several  feed-ports  designed  into  the  antenna 
system,  where  the  weighting  is  adjusted  according  to  some  specified  control 
algorithm.  For  communications  systems  where  the  antenna  field  of  view  (FOV) 
may  contain  both  interference  and  signal  sources,  adaptive  antennas  can  be 
utilized  to  improve  the  signal-to- interference  ratio.  This  is  done  by  auto¬ 
nomously  changing  the  pattern  so  that  relative  pattern  maxima  are  directed 
toward  the  desired  user  locations  and  pattern  minima  (nulls)  are  directed 
toward  the  interference  locations. 

The  basic  performance  of  such  an  adaptive  antenna  system  is  influenced 
somewhat  independently  by  two  sub-systems:  the  antenna  and  the  adaptive 
processor.  Choice  of  the  antenna  type  (multiple  beam  antenna  (MBA)  or  phased 
array  (PA))  and  configuration  (number  of  beams  or  array  element  positions)  is 
usually  determined  by  such  considerations  as  signal  interference  proximity 
(i.e.,  resolution),  number  of  interference  sources,  nulling  bandwidth,  number 
of  simultaneous  users  and  their  locations  over  the  FOV,  the  specific  FOV  over 
which  the  antenna  must  be  designed  to  operate,  etc.  Similarly,  the  design  of 
the  adaptive  processor,  which  is  used  to  control  the  weighting  at  each  ar.tenna 
feed-port,  depends  strongly  on  such  factors  as  how  the  user  signals  are  to  be 
distinguished  from  interference  signals,  the  time  required  for  the  processor 
to  "adapt"  to  a  steady  behavior,  the  bandwidth  of  operation,  etc.  For  these 
reasons,  the  specific  characteristics  of  existing  adaptive  nulling  systems 
are  as  different  as  there  are  applications  for  adaptive  nulling  techniques. 

To  date,  much  has  been  published  on  the  ideal  performance  characteristics 
of  adaptive  nulling  antenna  systems.  However,  little  has  been  published  on 
the  effects  of  hardware  component  imperfections  on  system  performance.  To 
characterize  these  effects  is  the  purpose  of  this  note.  We  present  methods 
for  categorizing,  analyzing  and  measuring  the  effects  of  a  few  of  the  key 
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components  which  are  common  to  the  majority  of  adaptive  systems  implemented 
to  date.  Some  examples  we  will  discuss  are  those  effects  due  to  antenna 
type,  antenna  tolerance  errors,  channel  tracking  errors,  weighting  circuits, 
hard- limiters,  correlators  and  base-band  components.  We  analyze  these  from 
the  viewpoint  of  the  degree  of  cancellation  achievable  by  the  adaptive  pro¬ 
cessor  and  also  their  impact  on  loop  dynamics  (where  appropriate).  We  will 
also  identify  several  techniques  we  have  found  useful  in  evaluating  these 
component  effects.  One  such  technique  has  already  been  published  in  Lincoln 
Laboratory  Technical  Note  1978-1^  and  we  will  utilize  this  material  where 
appropriate.  Others  will  be  developed  in  the  forthcoming  sections. 

The  characterization  of  the  performance  of  the  antenna  sub-system  has 
been  discussed  in  some  detail  in  Ref.  1;  hence,  the  major  thrust  of  this  note 
will  be  to  concentrate  on  those  factors  associated  with  the  adaptive  processor 
sub-system.  However,  a  brief  summary  of  antenna  dispersion  (i.e.,  frequency 
variations  introduced  by  the  antenna)  and  some  updated  estimates  on  the 
effect  of  antenna  type  (MBA  vs  PA)  on  cancellation  will  be  included  here.  In 
order  to  illustrate  the  various  component  effects,  measured  performance  results 
obtained  from  an  experimental  Applebaum-Howell  type  adaptive  processor  will 
be  presented. 
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1 1 .  BACKGROUND 


In  this  section  we  present  a  brief  summary  of  the  basic  equations  charac¬ 
terizing  the  performance  of  a  wideband  Applebaum- Howells  type  adaptive  nulling 
processor.  A  simplified  diagram  of  the  processor  is  illustrated  in  Fig.  1. 

The  narrowband  analysis  pertaining  to  this  type  of  processor  has  been  discussed 

(21 

in  some  detail  in  a  review  article  by  Gabriel  .  In  our  model,  we  include  the 
frequency  response  of  the  various  components  (e.g.,  the  antenna,  RF  front-ends, 
IF  filters,  etc.)  in  order  that  a  wideband  characterization  of  the  processor 
can  be  obtained.  We  also  include  in  our  analysis  the  effects  of  a  frequency- 
dependent  feedthrough  component  inherent  in  the  weight  circuitry.  We  adopt 
the  following  notation: 

n(t)  =  complex  envelope  of  a  single  interference  noise 
waveform  incident  on  the  aperture. 

Sn(u>)  =  power  spectral  density  of  the  noise,  n(t). 

Ak(a>)  =  frequency  response  function  of  the  kth  antenna 
element  (including  time  delay  in  the  case  of  a 
Phased  Array) . 

(w)  *  channel  frequency  response  function  charac¬ 

terizing  the  kth  channel  RF/IF  mixer,  amplifiers, 
and  filters,  over  the  nulling  bandwidth. 

E.  (u)  *  overall  frequency  response  of  kth  channel. 

Ek(*>)  -  Ak(«)  •  Hk(u.). 

v^(t)  =  complex  envelope  of  signal  voltage  at  the  input 

to  the  kth  weight  circuit. 

w* (t)  =  complex,  frequency  independent  weight,  where 

denotes  complex  conjugate. 

y£(w)  =  frequency-dependent  feedthrough  path  present  in 

the  weight  circuitry. 

v  (t)  =  complex  envelope  of  the  processor  output  voltage. 

X^(t)  =  complex  envelope  of  the  correlator  mixer  output. 

tq  =  time  constant  of  first  order  filter  following  the 
correlation  mixer. 


-  INCIDENT  INTERFERENCE 
WAVEFORM 


Fig.  1.  kth  channel  of  an  N-channel,  wideband  Applebaum-Howells 
type  adaptive  nulling  processor. 
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vk 


feedback  loop  gain. 

complex  beam  steering  voltage  specifying  the 
quiescent  (no  interference  sources)  weight 
setting  and  thereby  the  quiescent  radiation 
pattern. 

cross  correlation  between  vv(t)  and  v  (t). 


cross  power  spectral  density  of  v,  (t)  and  v  (t) 

K  Q 


wo  =  carrier  radian  frequency  of  input  waveform. 


BWc  =  closed- loop  bandwidth  in  hertz. 


BW  =  nulling  bandwidth  (as  determined  by  H,  (w)) 
in  hertz. 


R. 

1 


C 


(I/N) 


b,a 


total  power  of  n(t)  in  the  nulling  bandwidth 

Rj  =  17  f  Sn{to)  du  • 

cancellation:  C  *  (I/N)a/(I/N)b. 

ratio  of  interference  noise  level  to  thermal 
noise  power  before  adaption  (b)  or  after 
adaption  (a) . 


Complex  notation  is  used  because  most  commonly  employed  weighting  techniques 
use  the  inphase  (I)  and  quadrature  (Q)  outputs  of  the  correlation  mixer  to 
drive  separate  I  and  Q  weight  circuits  in  each  channel.  These  are  modelled 
as  a  single,  complex  weight  circuit  with  a  frequency  response  which  is 
given  by 

“k  +  Yk(u)  =  ak  ~  j®k  *  (Ykl('a>)  ~  jYk2(“))  Cl) 

where  and  6^  denote  the  I  and  Q  frequency  independent  weights  and  Ykl(w) 
and  Yj^C1*1)  denote  the  I  and  Q  frequency  dependent  weight  feedthrough  paths 
respectively,  and  6^  are  assumed  directly  controllable  over  the  range  -1  to 
+1  by  means  of  drive  voltages  derived  from  the  I  and  Q  outputs  of  the  correlator 
circuitry. 


5 


Using  the  above  notation,  the  processor  output,  vQ(t)  can  be  expressed 
as  a  weighted  sum  of  the  N  waveforms,  v^ft),  k  =  1,  ....  N: 

N  N  /*x> 

V  (t)  =  £  w*  v  (t)  J  T*  (-t)  v  (t-x)  dx  (2) 

q=l  4  4  q=l  -as  4  4 

where  we  have  defined  the  impulse  response  of  y(o»)  to  be  y(“>)  ■*-*•  r(t),  where 
the  arrow  denotes  a  Fourier  transform  pair.  Note  that  y(w)  r(t) 
implies  that  y*(“)  r*(-t),  leading  to  the  presence  of  r*(-x)  in  the 

convolution  integral  in  Eq.  (2).  By  considering  the  basic  feedback  rela¬ 
tions  defined  at  the  sum  junction  in  the  baseband  branch  of  the  loop,  it  can 
f  2  j 

be  shown  that  each  complex  weight,  w^  satisfies  the  first-order  differ¬ 
ential  equation, 

Vk  +  wk  +  uXk  =  Vk  ,  k  =  1,  ....  N  (3) 

where  X^(t)  is  the  output  of  the  kth  complex  correlation  mixer  and  the  dot 

notation  indicates  a  time  derivative.  The  complex  correlation  mixer  consists 
of  separate  I  and  Q  mixers  followed  by  low  pass  filters  to  select  the  lower 
side  band.  In  complex  notation  the  mixer  output  is  given  by  X. (t)  =  v,  (t)  v*(t). 
Note  that  X^(t)  is  a  stochastic  waveform.  It  is  conventional  to  design  the 
feedback  loops  so  that  the  closed  loop  bandwidth  is  considerably  less  than  the 
bandwidth  of  n(t)  after  filtering  by  the  bandpass  filters  in  the  RF/IF  front 
end.  In  this  case,  Eq.  (3)  acts  to  yield  an  averaging  or  integrating  effect 
on  xk(t).  the  input  process  is  stationary  over  the  loop  adaption  time, 
then  the  "short  time  average"  of  wk(t)  (averaged  over  the  time  scale,  l/BWc, 
where  BWc  is  the  closed-loop  bandwidth  in  hertz)  can  be  estimated  by  averaging 

Eq.  (3)  over  the  ensemble  of  sample  functions  of  n(t).  We  obtain 

Vk  +  wk  +  uE{XkJ  *  vk  ,  k  =  1,  . . . ,  N  (4) 

where,  for  convenience,  we  now  denote  by  the  same  variable,  w^.  It 

remains  to  determine  E(X^).  Using  Eq.  (2)  we  obtain 
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E{Xk)  =  l  w  E{v  Ct)v*(t)} 

nal  ^  ^ 


+  I  /”  r  (-t)  E{vk(t)v*(t-r)>  dr  (5) 

q=l  .00  H  K  M 

Consider  the  case  where  a  single  interference  source  is  incident  on  the  antenna. 
The  expected  value  in  Eq.  (5)  is  readily  expressed  in  terms  of  the  power 

spectral  density,  S  (u>)  of  the  interference  source  input  process,  n(t).  First, 

•  n  f  3") 

define  ^(t)  cross-correlation  function  between  vk(t)  and  v  (t) v  1 , 

and  denote  the  cross  power  spectral  density  to  be  S,  (w)  +-*  R  (t) .  Then 

k,q  K,q 

Eq.  (5)  can  be  written  as 

ElV  -  I,  ",  Rk,q(°>  *  !  r  r;(-0  <*’  <6> 

q=l  ^  q=l  -»  M 

The  integral  in  (6)  can  be  expressed  in  terms  of  S,  (a>)  by  replacing  r*(-x) 

K,q  q 

and  R,  (t)  by  their  Fourier  transforms. 

K ,  q 

r  rq(-0  S.qW  dT  ,  i -  r  r„M  st  ,(.)  d.  (7) 


2,  q  k,q 


Using  Eq.  (7),  Eq.  (6)  can  be  rewritten  in  the  form 

E{v  =  i  h  r  \ +  Tq(«))  sk  (u.)  du 

(Jr  |  —oo  1  ‘  y  ^ 

The  cross-spectral  density  S,  (<u)  is  readily  expressed  in  terms  of  S  (w) 

*  f  4  n 

using  the  results  of  Ref.  3  (Chapter  9): 


Sk,q(w)  =  V“)  Aq(w)  Hk(u)  Sn(w)  (9 

Finally,  we  assume  that  S^Coj)  is  white  over  the  nulling  band  and  define 
Sn(w)  =  Rj/BW  where  we  assume  each  Hk(w)  has  a  nominal  bandwidth  of  BW  Hertz, 
and  Rj  denotes  the  average  power  of  the  incident  noise  process  over  the 
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nulling  band  arising  from  a  single  interference  source.  In  this  case  Eq.  (8) 
can  be  written  in  the  simplified  form 


E{V 


N 

Ri  l  <(«  +  Y  (w))  Ak(u)  A* (w)  Hk(w)  H*(w)> 

J  n  =  ]  4  H  4  M 


(10) 


where  the  bracket  notation  <•>  denotes  an  average  over  the  nulling  bandwidth 

.)  dw  (11) 


u)  ♦itBW 
o 


,>  .  _L  f  , 

~  2ttBW  J  1 


co  -irBW 
o 


Using  Eq.  (10)  in  Eq.  (4),  the  time  evolution  of  the  N  weights  w^(t),  ....  w^,(t) 
can  be  expressed  in  vector  form  as 


t^w  +  w  +  y  £  Y?  •  (w  +  y)>  =  V  (12) 

where  the  notation  "t"  denotes  the  complex  conjugate  transpose  operation,  the 
underbar  denotes  a  column  vector  of  the  N  appropriate  scalar  variables,  and 
Ek(oo)  =  Hk(w)  Ak(io)  is  the  overall  channel  frequency  response. 

Equation  (12)  governs  the  time-evolution  of  the  complex  weight  w  for  a 
white  noise  process  as  an  input  to  the  processor  resulting  from  a  single 
interference  source,  and  with  frequency  dependent  channel  characteristics 
embodied  in  £(u) .  For  ideal  weight  circuits,  y(w)  =  0  and  Eq.  (12)  reduces  to 


t^w  +  w  +  yR  •  w  =  V  (13) 

where  the  matrix  R  =  as  defined  above.*  Equation  (13)  is  the  conven¬ 

tional  form  of  the  equations  characterizing  the  time-evolution  of  w  that  is 
discussed  in  Ref.  2  for  narrow-band  channels  and  in  Ref.  1  for  wide-band 
channels.  The  effect  of  y_  on  the  performance  of  the  processor  is  one  of  the 
component  error  effects  to  be  considered  in  the  following  sections. 


*A  more  general  definition  of  R  is 

R  =  E(v(t)  v(t)+)  (14) 

The  definition  used  in  Eq.  (13)  is  a  special  case  which  is  particularly  useful 
in  testing  or  characterizing  nulling  performance  in  the  frequency  domain. 
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We  conclude  this  section  by  defining  some  parameters  that  are  useful  for 
characterizing  the  performance  of  a  processor.  The  expected  value  of  power 
output  of  the  processor  is  given  by 

po  ■  WIV1'!2’  (15) 

Strictly  speaking,  this  is  an  ensemble  average  defined  over  the  ensemble  of 
sample  functions  of  the  input  noise  process,  n(t).  Practically  speaking,  it 
is  essentially  the  same  as  a  time  average  of  the  power  as  long  as  the  average 
is  taken  over  an  interval  of  time  greater  than  the  inverse  nulling  bandwidth 
(in  hertz).  Using  the  same  procedure  used  above,  Pq  can  be  expressed  as  a 
frequency  average 

PQ  =  <(w  +  1(W))+  •  Rj  E(u>)  E+(w)  •  (w  +  Y (<*>))>  (16) 

w  in  Eq.  (16)  is  frequency  independent.  If  f:(w)  and  y(id)  vary  significantly 
over  the  nulling  bandwidth,  the  average  output  power  will  necessarily  be 
greater  than  the  case  where  y_  and  E  are  frequency  independent.  To  assess 
these  effects,  it  is  convenient  to  define  the  cancellation,  C,  by 

C  =  (I/N)a/(I/N)b  (17) 

where  (I/N)  denotes  the  ratio  of  average  interference  output  power  (i.e., 

that  contribution  of  Pq  due  to  n(t))  to  the  average  thermal  noise  output 

power.  The  subscripts  "a"  and  "b"  denote  "after  adaption"  and  "before 

adaption,"  respectively.  C I/N) ^  is  a  function  of  the  steering  vector,  V, 

which  is  the  weight  setting  in  the  absence  of  interference  (or  quiescent 

weight  setting).  Following  Ref.  2,  it  is  conventional  to  normalize  the 

incident  interference  power,  R^ ,  to  receiver  thermal  noise  power.  When  this 

is  done,  the  cross- spectral  density,  S,  (u)  of  the  thermal  noise  is  an 

K,q 

identity  matrix.  This  assumes  the  noise  in  each  channel  to  be  uncorrelated 
with  the  other  channels  and  of  the  same  power  level  as  the  other  channels. 

Then,  using  Eq.  (16),  the  expression  for  (I/N)  ,  reduces  to 

a  OX  D 
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/ 


(I/N) 


a  or 


b 


<(w  ♦  y)+  •  R.  E  Ef 
^(w  +  y ) '  • 


•  (w  +  y)> 

(w  + 


(18) 


where  we  get  the  "after  adaption"  or  "before  adaption"  value  by  inserting 
the  appropriate  weight  vectors.  When  y  =  0  this  reduces  to  the  more  familiar 
form. 


(I/N) | 
Y=0 


t  n 

w  •  R  •  w 


t 

w  •  w 


(19) 


where,  as  before,  R  =  (E.  • 

The  above  definitions  are  the  basis  for  a  systematic  approach  to 
classifying  and  quantizing  the  effects  of  component  imperfections  on  the 
overall  performance  of  an  adaptive  antenna  system.  For  example,  antenna 
dispersion  is  accounted  for  by  A^(u).  Channel  tracking  errors,  or  the 
mismatches  in  the  channel  frequency  responses,  can  be  incorporated  in  the 
H^(u)  and  so  on.  In  the  following  sections,  these  error  effects  are  quanti¬ 
fied  in  terms  of  their  impact  on  the  achievable  cancellation,  C.  Cancellation 
is  not  entirely  satisfactory  for  use  as  a  figure  of  merit  because  it  depends 
not  only  on  the  error  effects,  but  also  on  the  choice  of  steering  vector,  V, 
power  level,  and  processor  threshold  level.  Quoting  a  cancellation  number 
alone  can  therefore  be  misleading  unless  one  specifies  the  other  conditions 
for  which  it  is  applicable.  We  have  generally  used  one  of  two  standard 
choices  of  beam  steering  vectors  in  the  following  discussions.  Unless 
otherwise  stated,  we  have  used  an  "earth  coverage"  quiescent  radiation 
pattern,  which  is  a  uniform  gain  beam  18°  wide.  The  corresponding  V  for  a 
PA  is  V  =  col(l,  0,  0,  ....).  The  other  standard  choice  is  a  V  that  points 
a  maximum  gain  beam  in  some  specified  direction. 


Finally,  before  proceeding  to  a  discussion  of  the  results  summarized  in 
Section  III,  some  basic  assumptions  inherent  in  these  results  should  be 
discussed.  First,  we  assume  weight  circuits  whose  phase  and  amplitude  are 
ideally  frequency  independent.  This  precludes  the  use  of  tapped  delay  line 
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or  other  frequency  dependent  weight  designs.  Second,  we  constrain  the 
processor  to  allocate  only  a  single  degree  of  freedom  per  interference 
source.  When  a  single  interference  source  is  incident  on  an  adaptive 
antenna,  the  adaptive  processor  will  generally  allocate  as  many  degrees  of 
freedom  as  are  needed  to  form  a  null  over  the  total  nulling  bandwidth. 

However,  this  situation  is  undesirable  for  several  reasons,  perhaps  the  most 
important  of  which  is  that  using  more  than  a  single  degree  of  freedom  to 
form  a  wideband  null  on  the  source  leads  to  a  spatially  broad  null.  (Note: 
for  an  array,  the  phase  at  each  element  varies  as  w  sine  --  or  simply,  as  w8 
for  small  0.  Hence,  the  processor  output  is  a  function  of  a> 6,  so  that 
fixing  0  and  varying  w  gives  the  same  result  as  fixing  u  and  varying  0.) 

Hence,  if  one  allocates  several  processor  degrees  of  freedom  to  null  a 
single  source  over,  say,  a  20%  frequency  band,  the  angular  width  of  the 
resulting  spatial  null  will  be  20%  of  the  interference  source  off-axis 
angle.  This  significantly  limits  the  off-axis  angular  resolution  achievable 
with  the  antenna.  Another  reason  for  constraining  the  processor  to  one 
degree  of  freedom  per  source  is  simply  that  the  systems  we  have  considered 
have  a  limited  number  of  degrees  of  freedom  and  we  want  to  maximize  the 
number  of  sources  they  are  capable  of  nulling. 

III.  CLASSIFICATION  OF  COMPONENT  ERROR  CHARACTERISTICS  AND  THEIR  EFFECTS 

A  block  diagram  illustrating  the  basic  sections  of  an  adaptive  antenna 
feedback  nulling  system  is  shown  in  Fig.  2.  We  use  a  feedback  nulling 
processor  to  develop  the  effects  of  the  various  components  on  performance. 

Most  results,  however,  are  readily  applicable  to  other  types  of  processors 
such  as  open-loop  and  feed-forward  processors.  The  feedback  mechanism,  however, 
does  compensate  for  certain  frequency  independent  errors  which  would  affect  the 
cancellation  performance  of  open-loop  systems.  Parameters  characterizing 
each  sub-section  of  the  processor  are  listed  in  each  block.  The  effects  of 
parameters  associated  with  each  sub-section  on  the  cancellation  achievable 
from  the  processor  will  be  developed  in  detail  in  Sec.  IV.  For  convenience 
these  results  are  summarized  in  Table  1.  A  brief  discussion  of  Table  1  will 
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TABLE  1 


FACTORS  AFFECTING  CANCELLATION 
EFFECT  CANCELLATION  DEPENDENCE 


1  .  Antenna  Dispersion 

a)  N-Element  Array 

b)  MBA  -  Time  Delayed  Beams 

2.  Correlated  Channel  Tracking  Errors 

(A  =  RMS  Path  Length  Mismatch 
consisting  of  Antenna  Placement 
Errors  Normal  to  the  Plane 
Containing  the  Array  or  Cable/ 
Channel  Path  Length  Mismatch) 

2 

3.  Uncorrelated  Channel  Tracking  Error  (o  ) 

a)  Earth  Coverage 

b)  Beam  to  User 

4.  Weight  Circuits 

2 

a)  Feedthrough  (Isolation  |y  |  ) 

°  2 

Uncorrelated  Feedthrough  (o  ) 
Correlated  Feedthrough  (delay  t) 

2 

b)  Weight  Quantization  (o  ) 

c)  I/Q  Imbalance  (Frequency  Independent) 

5 .  Errors  in  Feedback  Paths 

a)  Frequency  Independent 

b)  Frequency  Dependent 

c)  Weight  Noise 

6.  Hard-Limiting  in  Feedback  Loops 
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C  =  sin  6 *FBW^2 

C  =  K2(FBW)2 


A 


2 


C  —  °MAX 


C  -  °MAX  /N 


C  <  N | yq | 2  a2 

c  <  N2|T  I2 

C  <  a2 

No  effect. 


Equivalent  to  Beam  Steering 
Vector  Error  (Loop  Stability 
may  be  affected) 

Weak  Dependence 

Spectrum  Modulation  of 
Communications  Signals 

Small  Eigenvalue  Suppression 
Relative  to  Non-Limited  System 

s2  -►  s2//^ 


1  k 


be  presented  in  the  remainder  of  this  section;  the  reader  not  interested  in 
the  development  of  SectionlV  leading  up  to  these  results  may  readily  skip  to 
Section  V. 

Recognizing  that  the  results  of  Table  I  apply  with  the  assumptions 
mentioned  above  and  in  Section  II,  the  cancellation  dependence  on  each 
sub-section  of  the  processor  can  now  be  considered. 

1 .  Antenna  Dispersion 

The  effects  of  antenna  frequency  dependence  (i.e.,  antenna 
dispersion)  on  the  cancellation  achievable  with  a  particular  antenna  type 
have  been  studied  in  some  detail  in  Ref.  1.  Basically,  antenna  dispersion 
results  from  the  differing  frequency  responses  of  each  element  of  the  adaptive 
antenna.  Generally  each  antenna  output  port  has  a  frequency  response  function 
that  depends  on  the  interference  source  location.  The  nature  of  this  fre¬ 
quency  dependence  is  fundamentally  different  for  an  array  antenna  when  compared 
to  a  multiple  beam  antenna  employing  time  delayed  beams,  and  this  difference 
clearly  shows  up  in  the  dependence  of  cancellation  on  aperture  size.  Although 
the  cancellation  is  a  function  of  fractional  bandwidth  (FBW  =  BW/fQ, 

f  =  center  frequency,  BW  «  bandwidth  in  Hertz)  and  degrades  according  to 
°  2 

(FBW)  for  both  antenna  types,  the  cancellation  of  an  MBA  becomes  independent 

of  its  diameter  for  very  large  apertures  whereas  the  performance  for  the  PA 

2 

degrades  as  (D/X)  .  Therefore,  MBA's  become  inherently  more  broad  band  at 
higher  frequencies  when  larger  apertures  might  be  used.  A  theoretical 
development  with  computer  simulations  verifying  this  result  will  be 
presented  in  Sec.  IV. 

2.  Correlated  Tracking  Errors 

Assembly  of  the  antenna  system  usually  results  in  errors  in 
positioning  the  antenna  elements  and  matching  the  cable  lengths.  Generally, 
for  a  satellite  at  geosynchronous  orbit,  only  the  placement  errors  that  are 
perpendicular  to  the  plane  containing  the  array  are  important.  These  can  be 
described  statistically  as  an  rms  error,  A,  in  the  overall  path  length  of 
each  channel.  The  important  conclusion  from  the  result  of  Table  I  is  that 


matching  of  these  path  lengths  is  dependent  only  on  the  bandwidth  and  not  the 
RF  or  IF  frequencies.  The  cancellation  deteriorates  approximately  as  y(irBWA/Vp) 
where  v  is  the  phase  velocity  in  the  propagation  medium.  These  tolerances  on 
path  length  matching  must  be  maintained  following  the  antenna  through  each 


channel  down  to  the  sum  junction,  although  Vp 
propagation  medium. 


will  vary  dependent  on  the 


Uncorrelated  Channel  Tracking  Errors 


Implementation  of  the  N-channels  of  an  adaptive  processor 

usually  results  in  differences  in  the  frequency  response  functions  (i.e., 

tracking  errors)  of  each  channel.  These  tracking  errors  can  usually  be 

2  2  2 

characterized  by  an  rms  amplitude  and  phase  error,  o  =  °A  +  ,  which 

tenJs  to  be  uncorrelated  from  channel  to  channel.  The  tracking  errors  are 


generally  dominated  by  the  filters  in  each  channel  which  set  the  processor 
nulling  bandwidth  and  by  the  mixers  which  may  be  used  to  convert  each  channel 
to  an  IF  frequency  for  weighting.  The  effects  of  such  uncorrelated  tracking 
errors  have  been  discussed  in  some  detail  in  Ref.  1,  and  the  cancellation 
limits  for  two  particular  modes  of  operation  listed  in  Table  I  are  developed 
in  that  reference.  The  most  notable  of  those  results  is  that  in  the  case  of 


an  earth  coverage  quiescent  pattern  obtained  from  a  PA  antenna,  best  cancel¬ 
lation  results  are  obtained  by  selecting  the  "best"  channel  (i.e.,  the  one 
with  the  lowest  measured  rms  tracking  error)  as  the  reference  channel. 


Weight  Circuits 


The  usual  approach  to  adaptive  processing  is  to  implement 

frequency  independent  weights  --  i.e.,  ones  for  which  amplitude  and  phase 

hold  the  same  values  over  the  entire  nulling  h.nd.  In  reality  such  a  weight 

is  achievable  only  within  certain  limits.  We  have  found  that  a  good  model 

applicable  to  many  realistic  weight  circuit  designs  is  an  ideal,  frequency- 

independent  weight  shunted  by  a  non-weighted  feedthrough  path  (see  Fig.  1). 

2 

The  amount  of  feedthrough  is  characterized  by  an  isolation  parameter,  |ydI  , 
and  a  time  delay,  t,  relative  to  the  weighted  path. 
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This  model  can  accommodate  two  basically  different  types  of  feedthrough  -- 
correlated  and  uncorrelated.  In  the  latter  case  the  k**1  channel  feedthrough 
path  is  characterized  as  a  random  function  of  frequency,  y^  (<*>),  with  mean,  Yok» 
and  standard  deviation,  The  Y^C^)  are  assumed  to  be  uncorrelated  from 

channel-to-channel  in  the  sense  that 

J"  <(Yk(»)  -  Yok)  (Yj  (»)  -  Yoj)*>  o  =  0  (20) 

For  uncorrelated  feedthrough  we  assume  to  be  0. 

To  model  correlated  feedthrough  we  have  a  frequency- independent  isolation 

parameter,  y  in  series  with  a  time  delay  in  which  case  the  frequency  response 
t.h  ^  ^ 

of  the  k  1  feedthrough  path  is 

=  Yok  6  (21) 

The  y (uiD  are  obviously  correlated  from  channel  to  channel  with  the  case  of 
equal  time  delays  yielding  maximum  correlation.  The  correlated  feedthrough 
model  is  applicable  to  several  real  weight  circuit  designs  and  is  considered 
to  be  a  more  useful  model  than  uncorrelated  feedthrough. 

Notice  from  Table  I  that  the  cancellation  effects  of  feedthrough  are 

functionally  similar  to  the  effects  of  channel  tracking  errors  except 

2 

that  cancellation  performance  degrades  as  N  where  N  is  the  number  of 
channels.  When  the  processor  is  assembled,  two-channel  measurements  are 
the  simplest  to  make;  this  is  usually  the  manner  in  which  the  performance 

is  specified.  If  the  two-channel  measurement  yields  a  measured  cancellation, 

2 

C2,  then  the  system  performance  will  degrade  by  (N/2)  relative  to  C 2  in 
an  N-channel  system.  This  corresponds  to  a  12-dB  degradation  for  an 
8-channel  processor.  A  more  thorough  examination  of  feedthrough,  and 
techniques  for  measuring  this  effect,  is  presented  in  Sections  IV  and  V. 

Other  types  of  errors  which  might  occur  in  the  weighting  circuits  are 
I/Q  imbalance  and  weight-quantization,  if  digital  control  of  the  weights  is 
used.  I/Q  imbalance  results  from  an  I/Q-type  weight  where  the  phases  of  the 
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I  and  Q  sections  of  the  weight  are  not  in  true  quadrature  (for  example,  an 
imperfect  hybrid  is  used).  Normally  the  deviation  from  quadrature  is  on 
the  order  of  1-5°,  depending  on  the  IF  frequency.  However,  most  hybrids  are 
broadband  devices  so  that  this  deviation  from  quadrature  is  generally  constant 
over  the  nulling  band.  In  this  case,  it  can  be  shown  that  a  closed-loop  feed¬ 
back  algorithm  compensates  for  this  error  so  that  it  imposes  no  limitations  on 
cancellation.  If  the  I/Q  phase  difference  should  vary  with  frequency,  its 
effect  on  cancellation  will  be  the  same  as  for  a  weighted  channel -tracking 
error. 

Weight  quantization  usually  results  in  a  weight  setting  that  is  slightly 

different  from  optimum.  If  w  is  the  optimum  weight,  then  w^  ■  w  +  Aw  is 

the  applied  weight  as  a  result  of  quantization.  The  quantization  error,  Aw, 

is  similar  mathematically  to  an  rms  mismatch  between  channels  when  averaged 

2 

over  all  possible  weight  settings.  The  cancellation  is  then  limited  by  a  , 
which  we  defined  as  the  rms  quantization  error  (see  Eq.  (65)). 

5.  Feedback  Path  Errors 

The  basic  function  of  the  feedback  loops  in  an  adaptive  processor 
is  to  sense  the  N  x  N  correlation  matrix  of  the  N  signal  channels.  Denote  this 

correlation  matrix  by  R.  The  array  output  power  is  then  given  by  (without 

feedthrough) 

P  =  w^  •  R  •  w  (22) 

o  —  =  — 

where  w  is  a  column  vector  of  the  weights  and  "+"  denotes  the  complex  conjugate 
transpose  operation.  For  an  ideal  system,  then,  the  steady-state  adapted  weights 
are  given  by  (see  Eq.  (13),  with  tqw  =  0) 

w  =  [I  +  yR]"1  •  V  (23) 

where  y  is  the  effective  processor  gain  and  V  the  quiescent  beam  steering 
vector.  It  is  readily  shown  that  using  Eq.  (23)  in  (22)  results  in  a  minimum 
in  the  radiation  pattern  in  the  direction  of  all  interference  sources  while 


However, 


(4) 

yielding  the  best  rms  fit  to  the  quiescent  pattern  defined  by  V_. 
in  a  realistic  system,  due  to  the  inherent  errors  present  in  the  feedback  loops, 
a  slightly  erroneous  correlation  matrix,  R',  is  actually  sensed.  Two  basic 
questions  then  arise.  How  do  the  errors  in  R!  affect  the  achievable  cancel¬ 
lation  and  how  do  these  errors  affect  loop  stability  (i.e.,  the  transient 
response  during  adaption)?  A  stability  analysis  for  first-order  control  loops 
is  presented  in  Sec.  IV.  The  effects  on  cancellation  are  summarized  in  Table  I. 

The  effect  of  errors  in  R'  can  be  conveniently  separated  into  frequency 
independent  errors  and  frequency  dependent  errors.  Examples  of  the  former 
include  errors  in  quadrature  hybrids,  mixer  imbalance,  dc  offsets,  and  differ¬ 
ences  in  loop  gains,  whereas  examples  of  the  latter  would  include  loop  path 
length  mismatches,  and  loop  amplifier  and  filter  tracking  errors.  Frequency 
independent  errors  have  no  effect  on  the  capability  of  the  processor  to  null 
the  interference.  Their  main  effects  are  to  modify  the  equivalent  beam  steering 
voltages  and  to  modify  the  stability  margin  of  the  loops.  Frequency  dependent 
errors  in  the  feedback  loop  can  degrade  the  ability  of  the  processor  to  null 
the  interference,  but  such  effects  are  only  second  order  relative  to  channel 
tracking  errors.  c In  fact,  we  show  in  Section  IV  that  if  the  signal  channels 
track  perfectly,  then  frequency  dependent  errors  in  the  feedback  loops  are 
totally  compensated  for  by  the  feedback  mechanism. 

A  third  source  of  error  can  be  attributed  to  the  feedback  circuitry 
weight  noise,  either  generated  independently  in  each  loop  or  by  the  modulation 
induced  on  the  weights  as  a  result  of  the  interference. ^  The  main  effect  of 
this  "weight  jitter"  is  to  modulate  the  signal  waveforms  passing  through  the 
weight  circuits  making  the  processor  behave  as  a  time  varying  multi-path  com¬ 
munications  channel.  This  can  affect  the  signal  detectability  and  induce  cross 
talk  between  signals  closely  spaced  in  frequency.  For  example,  when  a  c-w 
signal  is  incident  on  the  antenna,  weight  noise  generated  internal  to  the 
baseband  circuits  modulates  this  signal  so  that  the  power  spectrum  of  the 
output  broadens,  resulting  in  a  higher  total  power  level  at  the  output.  This 
effect  can  be  used  as  a  diagnostic  tool  to  measure  weight  noise  arising  within 
the  feedback  loops. 


feedback  loops  has  led  to  the  conclusion  that  such  limiting  relaxes  circuit 
dynamic  range  requirements,  while  having  little  effect  on  the  system  cancella¬ 
tion  performance.  For  a  certain  class  of  feedback  loop  designs,  we  have  found 
this  result  to  be  valid  for  a  single  interference  source,  but  not  when  multiple 
interference  sources  are  present.  A  detailed  analysis  of  the  effects  of  hard- 
limited  nulling  processors  can  be  found  in  Ref.  7.  The  effect  of  hard-limiting 
on  cancellation  performance  (relative  to  a  non-limited  system  designed  for  the 
same  threshold  power)  is  summarized  in  Table  I  for  the  case  of  two  sources, 
and  which  generate  two  eigenvalues,  s^  and  s2,  of  the  correlation  matrix. 
Hard-limiting  in  the  feedback  loops  has  the  effect  of  increasing  the  sensitiv¬ 
ity  threshold  (or  equivalently  decreasing  the  feedback  loop  gain)  applied  to 
the  smaller  eigenvalue.  For  example,  if  the  larger  eigenvalue  is  40  dB  above 
the  design  threshold,  the  effective  threshold  for  the  smaller  eigenvalue  is 
raised  by  20  dB  and  the  smaller  eigenvalue  will  not  be  sensed  or  nulled  by  the 
processor  unless  it  is  20  dB  or  more  above  the  design  threshold.  The  result  is 
a  reduced  sensitivity  to  the  weaker  sources,  and  reduced  suppression  of  these 
sources  relative  to  a  non-hard- limited  system. 
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IV.  ANALYSIS  OF  COMPONENT  IMPERFECTIONS  ON  CANCELLATION  AND  LOOP  DYNAMICS 

The  previous  section  summarized  some  estimates  of  component  limitations 
on  system  cancellation  capabilities.  In  this  section  we  develop  these  esti¬ 
mates  in  detail  and  also  consider  the  effects  of  loop  component  errors  on 

loop  stability.  Uncorrelated  channel-tracking  errors  and  hard-limiting 

r i  7) 

effects  are  considered  elsewhere  *  and  will  not  be  discussed  further  in 
this  note.  The  remaining  component  effects  will  be  developed  in  this  Section. 

A.  Antenna  Dispersion  Effects 

In  this  Section  we  assume  the  antenna  and  processor  to  be  "ideal"  (in 
the  sense  that  no  errors  or  imperfections  are  present)  and  consider  the 
fundamental  limitations  on  cancellation  imposed  by  the  antenna  itself. 

The  achievable  cancellation  for  a  given  size  aperture  depends  strongly 

on  the  operating  frequency,  fractional  nulling  bandwidth  and  field-of-view. 

Furthermore,  cancellation  depends  strongly  on  the  choice  of  the  steering  vector 

that  governs  the  quiescent  radiation  pattern  (e.g.,  earth  coverage  or  directed- 

beam-to-user) .  A  rule  of  thumb  for  geosynchronous  satellite  antennas  is  that 

ideal  antenna  dispersion  will  limit  achievable  cancellation  to  about  40  dB 

(4) 

over  a  3%  bandwidth  with  a  D/A  of  10  ,  for  interference  sources  within  8°  of 

bores ight  using  one  degree  of  freedom  per  interference  source.  An  aperture 
diameter  of  10A  is  a  lower  limit  for  what  would  be  required  on  a  satellite- 
borne  antenna  operating  at  UHF.  Since  aperture  sizes  much  larger  than  10A 
will  be  desired  at  higher  frequencies  and  since  fractional  bandwidths  larger 
than  3%  might  be  employed  in  spread- spectrum  communications  systems,  it  is 
desirable  to  extrapolate  this  rule  of  thumb  to  larger  D/A's  and  to  larger 
bandwidths.  Since  the  cancellation  performance  of  multiple  beam  antennas  and 
phased  arrays  is  different,  we  consider  them  separately.  We  shall  show  that 
the  MBA  has  considerably  greater  bandwidth  potential  than  the  PA  as  D/\ 
increases . 
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1 .  Phased  Array 

The  dependence  of  cancellation  on  bandwidth  and  aperture  size  can  be 
determined  directly  from  the  eigenvalues  of  the  expected  value  of  the  cor¬ 
relation  matrix  of  the  signals  at  the  antenna  output  ports.  Denote  this 
matrix  by  R,  the  corresponding  eigenvalues  by  s^,  s2,  ....  s^  and  the  eigen¬ 
vectors  as  <e . ,  e?,  e^.  Then  the  cancellation,  C,  is  generally  given 

by^ 

C  =  K  S2/S1  (24) 

where  K  is  a  constant  depending  on  the  steering  vector  V  and  array  configuration 

+  "f 

and  we  assume  that  •  V  /  0  (the  case  e_2  *  X  =  ®  would  correspond  to  a  special 
case  unique  to  a  particular  choice  of  V  --  in  this  case  Cots^/s^).  It  remains 
now  tc  estimate  Sj  and  s2  as  functions  of  D/X  and  bandwidth.  In  Appendix  I  we 
show  that  for  the  phased  array, 

C  =  K  (D/X  sine  *  FBW)2  (25) 

where  K  is  a  constant  and  0  is  the  scan  angle  relative  to  array  broadside. 

This  is  the  functional  dependence  shown  in  Table  I.  For  example,  since  C 
2 

increases  as  (D/X)  ,  increasing  D  from  10X  to  100X  for  a  fixed  antenna  geometry 

degrades  the  achievable  nulling  performance  by  20  dB.'  For  D/X  =  10,  0  =  8° 

-4 

and  FBW  =  .032,  a  typical  number  for  C  is  10  and  we  can  express  (25)  in  the 
form 

C  =  5.2  x  10' 2  (D/X  sine  *  FBW)2  (26) 

Of  course,  better  cancellation  levels  might  be  obtained  by  using  up  more  than 
one  degree  of  freedom  per  interference  source,  but  this  results  in  a  broader 
angular  null  and  therefore  less  gain  for  users  in  close  proximity  to  the  inter¬ 
ference  source.  It  also  limits  the  maximum  number  of  interference  sources 
which  can  be  nulled  by  the  array.  Equation  (26)  indicates  that  for  a  frac¬ 
tional  bandwidth  of  2%,  null  depths  on  the  order  of  25  dB  are  realizable  for  an 
aperture  diameter  of  100X,  using  a  single  degree  of  freedom  per  interference 
source.  A  wider  fractional  bandwidth  of  5%  would  reduce  the  cancellation 
performance  to  less  than  -20  dB  for  the  same  size  aperture. 
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As  an  example  to  illustrate  the  validity  of  the  approximations  leading 
up  to  Eq.  (25),  consider  the  13-element  thinned  array  illustrated  in  Fig.  3 
having  the  normalized  coordinates  for  each  element  indicated  in  the  Figure. 

A  signal  source  is  positioned  at  angular  coordinates  6g  =  7.7°,  <J>s  =  0°  and 
an  interference  source  0.3°  away  at  0^  =  8.0°,  <(>.  =  0°,  which  is  equivalent 
to  an  elevation  angle  to  the  satellite  of  20°.  This  spacing  of  0.3°  between 
user  and  interference  corresponds  to  a  half-power  beamwidth/2  for  D/X  %  100. 
Figure  4  illustrates  the  behavior  of  the  first  four  eigenvalues  of  the  correla¬ 
tion  matrix  of  interference  signals  as  functions  of  D/X  assuming  a  fixed 
fractional  bandwidth  of  .046  and  white  noise  interference  spectrum  over  this 
band,  s^  corresponds  to  the  array  output  power  realized  with  a  maximum  gain 
beam  pointed  to  the  interference  source.  Therefore,  we  expect  Sj  to  be  very 

nearly  equal  to  R.  NG  (where  G  is  the  element  gain)  and  to  be  independent  of 
j  e  e 

D/X.  Observe  from  the  figure  that 


S1 

constant 

(27a) 

S2 

'V, 

(D/X  sin0 

*  FBW) 2 

(27b) 

S3 

(D/X  sin0 

*  FBW)4 

(27c) 

S4 

(D/X  sin0 

*  FBW)6 

(27d) 

The  behavior  of  s,  is  clearly  consistent  with  the  results  of  Appendix  I. 

To  check  the  validity  of  Eq.  (26)  for  cancellation,  we  compute  from  our 
example  the  cancellation  and  user  gain  loss  assuming  that  the  interference 
power  is  30  dB  above  the  processor  threshold.  This  threshold  level  of  the 
processor  is  indicated  in  Fig.  4.  Since  the  processor  does  not  respond  to  the 
eigenvalues  below  threshold,  only  eigenvalues  above  threshold  are  sensed  and 
nulled.  The  number  of  eigenvalues  above  threshold  is  the  number  of  degrees  of 
freedom  used  by  the  processor  to  form  the  null.  The  loss  in  user  gain  and 
cancellation  are  shown  in  Fig.  5  as  a  function  of  r\/X.  We  have  assumed  a 
steering  vector  that  directs  a  maximum  gain  beam  to  the  user  in  the  unadapted 
case.  The  user  power  is  negligible  in  its  impact  on  processor  performance. 
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Fig.  3.  13-element  thinned  planar  array  geometry  used  for  the 
simulations  of  Figs.  4  and  5. 
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Fig.  4.  Eigenvalues  Sj,  S2,  S3  and  sa  vs.  aperture  diameter, 
D/A  for  the  13-element  planar  array  of  Fig. 3.  The  interference 
source  is  positioned  at  0.  *  8°,  ij>.  =  0°.  FBW  =  .0455. 
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Fig.  5.  Cancellation  and  user  loss  in  directive  gain  vs.  DA 
when  a  user  source  is  positioned  0.3°  from  the  interference 
source  location  of  Fig.  4  for  the  13-element  planar  array  of 
Fig.  3.  A6  =  0.3°  corresponds  to  a  quarter  half-power  beam- 
width  (HPBW)  for  D/X  %  100. 
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Observe  that  for  D/X  £  15,  the  processor  senses  only  the  first  eigenvalue, 
therefore  only  one  degree  of  freedom  is  used  to  null  the  interference.  The 
resulting  cancellation  is  about  -30  dB.  For  D/X  -  25,  the  second  eigenvalue 
is  sensed  so  that  two  degrees  of  freedom  are  used  to  form  a  wide-band  null. 

The  resulting  cancellation  is  about  -25  dB,  4  dB  worse  than  predicted  by 
Eq.  (26)  (which  assumes  only  a  single  eigenvalue  sensed).  As  D/X  increases 
still  further,  the  cancellation  remains  fairly  constant  at  about  -28  dB,  even 
though  the  processor  is  now  using  more  degrees  of  freedom  to  form  a  null. 
However,  the  user  gain  degrades  rapidly  as  a  consequence  of  the  loss  in 
resolution  due  to  the  wide-band  null.  To  illustrate  this,  the  loss  in  gain  to 
the  user  which  would  be  obtained  for  FBW  =  0  is  also  plotted  in  the  figure. 

Up  to  7.5  dB  loss  in  directive  gain  in  the  direction  of  the  user  occurs  as  a 
consequence  of  using  additional  degrees  of  freedom  to  wide-band  the  null. 
Therefore,  the  achievable  S/J  (signal-to-interference  power  ratio)  improve¬ 
ment  averaged  over  the  band  is  only  about  16  dB  for  D/X  =  100! 

These  results  show  that  Eq.  (26)  is  valid  as  long  as  a  single  eigenvalue 
per  source  is  sensed.  They  also  illustrate  that  using  more  than  one  degree  of 
freedom  may  improve  the  cancellation,  but  it  may  not  improve  the  S/J  perform¬ 
ance  because  the  user  gain  may  suffer.  These  results  also  show  that  the  phased 
array  becomes  relatively  narrow-band  near  the  edge  of  the  FOV.  The  cancella¬ 
tion  versus  bandwidth  achievable  with  an  MBA  is  much  better  than  for  a  PA.  We 
now  consider  this  antenna  configuration. 

2 .  Multiple  Beam  Antennas 

The  general  characteristics  of  multiple  beam  antennas  for  use  in 
satellite-earth  coverage  applications  are  discussed  in  Ref.  (8).  The  salient 
features  of  MBA's  used  for  this  application  briefly  be  summarized. 

A  multiple-beam  antenna  located  on  a  geosynchronous  satellite  illuminates 
the  earth  with  a  fixed  position  set  of  beams,  the  composite  of  which  cover  the 
earth  FOV.  For  earth  coverage  illumination  (i.e.,  a  circular  FOV),  the  beams 
are  generally  positioned  hexagonal ly  and  usually  are  chosen  to  number  7,  19, 

37,  61,  ....  in  order  to  fully  illuminate  the  circular  FOV.  To  each  beam-set 
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covering  the  earth,  there  corresponds  a  given  size  aperture  consistent  with 
this  number  of  beams.  When  a  lens  or  reflector  antenna  system  is  used, 
physical  realization  of  the  beam-set  is  obtained  by  offsetting  an  array  of 
feeds  from  the  focal  point  of  the  lens  or  reflector,  transverse  to  the  focal 
axis.  Consequently,  each  beam  has  a  beam-width  and  sidelobe  levels  consistent 
with  full  aperture  illumination.  As  the  number  of  beams  increases  to  cover  a 
fixed  FOV,  the  aperture  size,  D/A,  must  increase  correspondingly.  Hence,  for 
a  large  number  of  beams,  there  is  significant  coupling  only  between  adjacent 
beams.  This  leads  to  a  correlation  matrix  R,  defined  at  the  feed  port  outputs, 
which  is  diagonally  dominant.  Furthermore,  an  interference  source  positioned 
at  a  specific  location  on  the  earth  is  characterized  only  by  that  set  beams 
in  close  proximity  to  the  source.  Consequently,  the  eigenvalues  of  the  feed- 
port  correlation  matrix  resulting  from  this  source  are  dominated  by  the  sub¬ 
matrix  generated  by  the  beams  surrounding  the  interference.  Otherwise  said, 
adding  more  beams  to  the  antenna  configuration  far  away  from  this  source 
should  not  affect  the  eigenvalues  of  R.  However,  as  discussed  above,  adding 
more  beams  corresponds  to  increasing  D/A  for  a  fixed  FOV.  Consequently,  we 
anticipate  that  the  cancellation  dependence  of  the  MBA  operating  over  a  fixed 
FOV  becomes  independent  of  D/A  (or  number  of  beams)  as  D/A  becomes  large. 

Thus,  the  cancellation  achievable  with  an  MBA  tends  to  be  dependent  on  the 
position  of  an  interference  source  within  the  surrounding  beam  set,  and  this 
positioning  is  predominantly  aperture  independent.  The  resulting  simulations 
will,  indeed,  verify  this  conjecture. 

Because  the  frequency  transfer  function,  A^(io),  characterizing  the  N 
beams  of  the  MBA  is  much  more  complicated  than  for  the  phased  array,  it  is 
difficult  to  obtain  analytical  results  predicting  the  dependence  of  cancel¬ 
lation  vs.  D/A  sinQ  and  FBW  for  the  MBA  similar  to  those  of  Appendix  I  for 
the  phased  array.  Consequently  we  resort  to  computer  simulations.  In  order 
to  generate  a  simulation  for  the  MBA  paralleling  the  analysis  of  the  phased 
array,  we  again  separate  the  jammer  and  user  positions  by  0.3°.  The  MBA 
performance  curves  corresponding  to  Figs.  4  and  5  for  the  phased  array  are 
illustrated  in  Figs.  6  and  7.  For  D/A  =  10,  7  beams  cover  the  FOV  and  were 


27 


EIGENVALUES  (dB)  -  RELATIVE  SCALE 


LOSS  IN  DIRECTIVE  GAIN  TO  USER 


0 

10 

20 

30 

40 

50 

60 


Fig.  7.  Cancellation  and  user  loss  in  directive  gain  vs.  D/x  for 
an  MBA  employing  time-delayed  beams  when  a  user  source  is  positioned 
0.3°  from  the  interference  source  location  of  Fig.  6.  A6  =  0.3° 
corresponds  to  (HPBW)/2  for  D/X  £  100. 
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used  in  the  simulation.  For  D/X  K  20,  the  simulation  employed  only  the  19  beams 
at  the  center  of  the  FOV  so  that  existing  computer  programs  limited  to  a  small 
number  of  beams  could  be  used.  The  user  and  interference  source  were  placed  at 
9g  =  1°,  4>s  =  0  and  =  1.3°,  4k  =  0  respectively.  Since  the  MBA  is  scan 
angle  dependent  only  relative  to  the  position  of  the  interference  source  amongst 
a  triad  of  adjacent  beams,  this  positioning  relatively  close  to  array  broad¬ 
side  does  not  compromise  the  simulation  results  concerning  bandwidth.  Inclus¬ 
ion  of  the  outer  beams  to  cover  the  FOV  completely  would  not  affect  the  results 
appreciably.  Observe  in  Fig.  6  that  all  the  eigenvalues  are  relatively  inde¬ 
pendent  of  aperture  size.  Generally,  for  fixed  input  power,  s.  should  increase 

2  1 
as  (D/X)  ,  since  Sj  represents  the  maximum  power  received  with  a  maximum 

directivity  beam  pointed  to  the  source,  and  this  maximum  gain  increases  as 
2 

(D/X)  .  However,  in  order  that  the  nulling  processor  used  in  the  simulation 
might  see  the  same  input  power  for  each  D/X,  we  have  adjusted  the  input  power 
for  each  D/X  so  that  Sj  remains  constant,  independent  of  D/X.  The  processor 
gain,  p,  is  then  chosen  for  each  case  so  that  ps^  =  30  dB.  The  simulations 
confirm  that  s2  is  also  independent  of  D/A.  Then  to  first  order,  the  cancel¬ 
lation,  C,  is  also  independent  of  D/A  because  C  «  s2/Sj  (Ell-  (24)).  Further¬ 
more,  we  note  that  Sj  is  relatively  independent  of  FBW  (s ^  corresponds  to  the 
output  power  when  a  maximum  directivity  beam  is  pointed  to  the  interference). 

This  beam  is  independent  of  frequency  for  the  small  percentage  bandwidths 

2 

(£  5%)  of  interest  here)  whereas  s2  increases  as  (FBW)  .  Hence, 

C  =  K  s2/Sj  =  K  (FBW)2  (28) 

For  an  MBA,  therefore,  C  is  predominantly  a  function  only  of  FBW.  This  is  the 
result  referred  to  in  Table  I.  The  MBA  cancellation  is  illustrated  in  Fig.  7, 
and  is  considerably  better  than  that  for  the  PA,  averaging  about  40  dB  over 
the  range  of  D/X  considered.  Furthermore,  the  cancellation  is  achieved  using 
only  a  single  degree  of  freedom.  This  is  evidenced  by  examining  the  directive 
gain  in  the  direction  of  the  user  assuming  both  FBW  =  0  and  FBW  =  4.55%,  which 
yields  similar  results  for  each  FBW.  For  D/X  =  100,  the  loss  in  user  gain  at 
one-half  power  beamwidth  from  the  interference  source  is  5  dB  rather  than  the 
predicted  3  dB  for  a  phased  array  when  FBW  =  0.  This  is  because  the  resolution 
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of  the  MBA  is  poorer  when  the  user  is  not  positioned  at  a  beam  center.  A 
maximum  gain  beam  synthesized  from  a  dyad  or  triad  of  beams  has  a  broader 
beamwidth  at  the  edge  of  the  19-beam  FOV.  Clearly,  however,  the  overall  S/J 
improvement  is  good  compared  to  the  phased  array.  For  example,  for  D/X  =  100, 
the  S/J  improvement  is  33  dB  for  the  MBA  compared  to  only  16  dB  for  the 
phased  array.  We  should  mention  again  that  the  results  of  the  simulation  are 
steering  vector  dependent,  so  that  deeper  cancellation  levels  might  be 
obtained  for  other  steering  vectors. 

One  final  note  which  should  be  mentioned  regarding  the  MBA  is  that  for 
employing  an  MBA  with  aperture  sizes  D/X  =  100,  and,  with  adjacent  beams 
spaced  one  HPBW  apart,  approximately  721  feeds  are  required  for  full  earth 
coverage.  Clearly,  employing  full  adaptive  control  on  each  feed  is  imprac¬ 
tical.  One  technique  for  simplifying  the  hardware  would  be  to  employ  a 
smaller  feed  cluster  which  could  be  scanned  to  cover  particular  areas  within 
the  FOV.  Another  technique  might  employ  a  switch  network  to  utilize  only 
those  beams  containing  signals  in  an  adaptive  processor  with  a  smaller  number 
of  channels. 

B.  Element  Positioning  and  Cable  Length  Mismatch 

An  array  antenna  in  space  may  experience  variations  in  its  element  posi¬ 
tions  or  electrical  path  lengths  as  a  result  of  the  extreme  temperature  vari¬ 
ations  to  which  the  outer  surface  of  a  satellite  is  subjected.  The  effect  of 
such  variations  on  cancellation  will  be  analyzed  in  this  section.  If  A  is 
the  rms  differential  displacement  of  the  elements  (expressed  in  feet),  then 
for  the  cancellation  to  be  a  level  C  over  a  bandwidth  BW  (in  MHz),  the 
tolerance  on  A  will  be  shown  to  be 

A[feet]  <  542.6  /c/BW  (29) 

Some  values  of  A  for  selected  bandwidths  and  cancellation  levels  are  tabulated 
in  Table  2.  It  should  be  emphasized  that  these  tolerances  on  A  assure  that 
the  antenna  element  position  errors  will  not  degrade  cancellation  below  the 
value  shown.  However,  they  do  not  assure  that  cancellation  to  these  levels 
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will  actually  be  achieved  because,  for  some  scenarios,  the  inherent  dispersion 
of  the  antenna  due  to  the  wavefront  angle  of  arrival  dominates  the  cancellation 
performance  so  that  in  those  cases  position  errors  are  not  important. 

TABLE  2 

REQUIRED  TOLERANCES  ON  RMS  DIFFERENTIAL  ELEMENT  DISPLACEMENT 
A  VS.  BANDWIDTH  AND  CANCELLATION  LEVEL 


BW (MHz) 

C(dB) 

A  (ft) 

A  (in) 

1 

-30 

17.2 

1 

-40 

5.4 

1 

-SO 

1.7 

10 

-30 

1.7 

10 

O 

l 

6.5 

10 

i 

cn 

O 

2.1 

100 

-30 

2.1 

100 

-40 

0.65 

100 

-SO 

0.21 

1000 

-30 

0.21 

1000 

1 

-u 

o 

0.07 

1000 

-50 

0.02 

Development 

Let  (xky^,0)  be  the  desired  x,y,z  coordinates  of  the  kth  element 
normalized  to  D/2,  where  D  is  the  smallest  diameter  enclosing  the  array.  Let 
the  matched  cable  length  extending  from  the  processor  of  each  element  be  £. . 
When  the  element  is  displaced,  and  the  cable  expands,  the  new  values  of  these 
variables  are  (x^  +  Ax^,  y^  +  Ay^,  Az^)  and  l  *  Afc^,  respectively.  Denote 
the  output  signal  at  the  kth  port  with  perfect  positioning  by  E^°.  Then,  for 
a  wavefront  of  amplitude  Eo  incident  from  angle  (0,<{>)  with  respect  to  the 
array  reference  frame,  E^°  is  given  by 
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(30) 


(«) 


j  kQ  ^  w(xjc  sin6  cos$  +  yk  sin0  sin$)  -jk  HW 

E  e  e 

o 


where,  as  before,  w  is  the  actual  angular  frequency,  is  the  angular  frequency 
at  band  center,  kQ  =  wq/c,  and  W  =  w/uq.  When  the  elements  are  displaced  and 
the  cables  expanded,  then  the  output  E^(uj)  becomes 


-jk  At.  W  +  jk  j  W[Ax,  sin©  cosi>  +  Ay,  sin0  simf> 

E/--»  OK  O  ^  K  K 

k(fc>)  =  Ej.  (“)  e  e 


+  Az^  COS0 ] 


(31) 


Note  that  errors  Ax^,  Ayfc  are  reduced  in  magnitude  by  sin0  (sin©  <  .15  within 
the  8.7°  FOV)  so  that  the  co-planar  positioning  of  the  elements  is  less  critical 
than  out  of  plane  position  or  cable  length.  Furthermore,  for  some  angles  of 
arrival,  the  dispersion  inherent  in  E^°(u)  will  dominate  all  positioning  errors, 
particularly  for  large  antennas.  However,  if  we  wish  to  specify  the  antenna 
tolerances  so  that  their  effect  on  cancellation  is  negligible  for  all  angles 
of  arrival,  then  the  case  of  sin0  =  0  (broadside  incidence)  is  a  worst-case 
condition. 


Define  M,  (oj)  =  E,  (w)  E*(ui).  Then  M  takes  the  form 
K,q  K  q  - 


-jk  (At,  -AH  )  W  +jk  -^(Az, -Az  )  W 
M  (a,)  =  M°v  (m)  e  0  k  <  e  °  2  k  W 

=M  =  K’q 


(32) 


where  we  have  used  the  small  angle  approximation  cos0  =  1  and  defined 

M.  (w)  =  E,°(u)  E  °*(u)).  Equation  (32)  indicates  that  the  dominant  errors 
=  K ,  q  K  q 

are  cable  expansion  and  deviations  in  out-of-plane  element  positioning. 
Define  the  net  position  error  for  the  kl  element  to  be 


Atk  -  AzkD/2 


Then,  assuming  broadside  incidence  ^Mk  ^(u>) 

_  -jk  W(A.  -A  ) 

1 2  _  J  o  v  k  q 


|  Eq  |  2) ,  we  obtain 
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(33) 


(34) 
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In  order  to  compute  the  effect  on  the  cancellation  level,  we  average  M(u>) 
over  the  frequency  band.  Define  R  =  (M(w.)^  ,  where  denotes  a  frequency 
average  over  the  bandwidth  of  interest:  Now  expand  R  in  the  form 

R  =  R  +  AR  (35) 

=  =o  = 

where  Rq  =  <Mq(u))^  and  AR  denotes  a  perturbation  of  Rq  due  to  the  net  position 
errors  of  each  element.  For  a  single  interference  source  incident  on  the  array, 
it  can  be  shown ^  that  the  cancellation  C  for  an  earth  coverage  beam-steering 
vector  is  bounded  by 


C  <  N(s2/Sj) 


(36) 


where  N  is  the  number  of  elements  and  s^  and  s^  are  the  two  largest  eigenvalues 
of  R.  (Normally  s^  and  s2  is  determined  by  the  channel  tracking  errors 

or  differences  in  frequency  response  among  channels,  which  in  this  case  result 
from  the  net  positioning  errors.)  To  determine  s^  and  s.,,  we  average  F.q.  (34) 
over  FBW  and  obtain 


R 


(A.  -A  ) 
k  q' 


(371 


To  express  (37)  in  the  form  of  Eq.  (35),  assume  that  kQ|A^-A  |  <<  1  and  also 


k  FBW 

that  — — ? —  |  A.  -  A  I 
2  1  k  q ' 


|  A,  -  A  |  <<  1 ,  where  BW  is  the  nulling  bandwidth  in 

C  K  C| 


hertz.  Then  (37)  can  be  written  in  the  form 

S  ■  lEol2  j1  'W2  *  jko(VV 


(38) 


It  is  reasonable  to  ass'ume  that  each  Ak  is  an  independent  random  variable  with 
zero  mean.  Let  A  be  the  standard  deviation  of  the  distribution  of  position 
errors,  then 


0  for  k=q 
2A^  for  Mq 


(39) 
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where  £{•}  denotes  the  expected  value  taken  over  the  distribution  of  position 
errors.  The  expected  value  of  R,  denoted  by  R,  is  therefore  given  by 


=k,q 


1 


_  i_  /rcBWx 
3  V  c  ) 


A 


(40 


where  6^  is  the  Kronceker  delta  function.  Equation  (40)  can  be  expressed  in 


the  form 


R  = 


VI  { 


(1-e)  1  1  +  el 


(41 


where  we  have  defined  1  =  col  [1 , 1 ,  .  .  . ,  1  ] ,  e  =  (“-)  A  >  "+"  is  the  complex 
conjugate  transpose,  and  I  is  the  identity  matrix.  The  eigenvalues  s^  and  s7 
of  R  can  readily  be  determined  to  be 


Sj  =-|Eor(l-c)  N 


(42 


S2  =  l£ol  e 


so  that,  using  Eq.  (36)  bounding  the  cancellation  C,  we  obtain 

2  „ 


C  < 


e  s  1  /ttBWV 

1-e  3 V  c  ) 


(43 


Expressing  BW  in  MHz  and  A  in  feet  then  leads  to  Eq.  (29)  specifying  the 
allowable  tolerances  on  A  as  a  function  of  the  desired  cancellation  and 
bandwidth.  Recall  that  (43)  is  based  on  a  statistical  average.  There  will 
be  some  variation  in  the  actual  cancellation  achieved. 


C.  Weighting  Circuits 

There  are  several  methods  for  implementing  amplitude  and  phase  control 
in  each  channel  of  an  adaptive  processor.  Perhaps  the  simplest  technique  is 
an  attenuator  and  phase  shifter  in  series  as  illustrated  in  Fig.  8a.  This 
topology  can  be  quite  useful  for  certain  applications.  Its  limitations  lie 
in  its  inability  to  turn  off  completely,  due  to  the  finite  dynamic  range  of 
real  attenuators,  but  perhaps  more  significantly,  it  is  difficult  to  couple 
this  type  of  weight  to  the  output  of  an  analog  adaptive  processor.  As  shown 
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(a) AMPLITUDE  AND  PHASE  WEIGHTING 


Fig.  8.  Two  common  techniques  used  to  obtain  amplitude  and  phase 
control  at  the  output  of  each  antenna  port:  (a)  Attenuator  in 
series  with  phase-shift,  (b)  Inphase  (I)  and  quadrature  (Q) 
weighting. 
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in  Section  III,  the  weight  control  output  of  the  processor  is  usually  avail¬ 
able  as  "inphase"  (I)  and  "quadrature"  (Q)  voltage  for  each  weight.  This 
output  must  be  converted  to  an  amplitude  and  phase  form  to  drive  the  attenu¬ 
ator  and  phase  shifter,  thus  complicating  the  circuitry.  To  avoid  this,  the 
weighting  is  usually  accomplished  in  terms  of  the  1  and  Q  components MiTectly, 
as  illustrated  in  Fig.  8b.  In  this  case,  the  I  and  Q  components  are  separated 
by  a  quadrature  hybrid.  The  I  and  Q  branches  can  then  be  driven  directly  by 
the  I  and  Q  components  of  the  correlator  output.  When  implemented  as  in  Fig. 

8b,  each  complex  weight  circuit,  depending  on  the  weight  setting,  can  have  as 
much  as  12  to  14  dB  of  insertion  loss  from  the  splitting  and  combining  losses 
of  the  hybrids.  A  lower  insertion  loss  topology  for  I  and  Q  weighting  is 
illustrated  in  Fig.  9.  In  this  scheme,  both  attenuation  and  phase  reversal  are 
accomplished  with  two  fewer  hybrids  than  in  Fig.  8b,  resulting  in  a  6  dB 
improvement  in  insertion  loss.  The  drive  voltage,  V^,  which  is  derived  from 
the  Q  control  voltage,  controls  the  pin-diode  resistance,  which  varies  the 
reflection  coefficients  F^  and  V ^  (which  are  nominally  equal)  over  the  range 
-1  to  +1.  This  allows  both  attenuation  and  phase  reversal.  A  direct  conse¬ 
quence  of  this  technique,  however,  is  the  presence  of  a  "feedthrough"  compo¬ 
nent  which  is  independent  of  control  voltage.  In  this  case  it  arises  in  the 
hybrid  parasitic  coupling.  In  practice,  a  small  component  of  feedthrough  is 
usually  present  in  any  weighting  circuit.  It  has  an  important  effect  on 
performance  and  it  is  of  interest  to  analyze  its  effect  on  cancellation. 

Development 

The  basic  equations  characterizing  the  effects  of  weight  feedthrough 
on  the  adapted  output  power  have  been  developed  in  Section  II.  The  composite 
weight,  including  feedthrough,  is  modeled  as  a  single,  complex  weight  circuit 
with  a  frequency  response  which  is  given  by  Eq.  (1).  That  part  of  the  weight 
which  is  directly  controllable,  denoted  by  w,  via  the  I  and  Q  correlator  out¬ 
puts  is  governed  by  Eq.  (12).  The  corresponding  output  power  accounting  for 
the  presence  of  feedthrough  is  given  by  Eq.  (15).  These  equations  are  generally 
difficult  to  solve  unless  one  resorts  to  numerical  simulations.  However,  if 
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Fig.  9.  Low  insertion  loss  weight  topology  for  obtaining  amplitude 
and  phase  control. 
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we  assume  an  array  antenna  geometry,  and  broadside  incidence  for  the  interference 
source,  then  the  questions  simplify  to  a  tractable  form.  With  these  assumptions, 
E^  E^  is  independent  of  frequency  and  Eqs.  (12)  and  (15)  reduce  to  the  simpler 
form 


T0W  +  [I  +  PR]  *  (w  +  <X>)  =  V  +  <y> 


(44) 


and 


?  =  ^  (w  +  <X>  +  AX)+  ‘  E  '  (x  +  <X>  +  ^x)  ^ 


(45) 


where,  as  before,  we  have  defined  R  =  R.  EE  and 

1 - 


y(w)  =  <y(w)>  +  Ay(w) 


(46) 


Note,  by  definition,  <&y>  =  0.  In  the  steady  state,  solutions  to  Eq.  (44) 
take  the  form 

w  +  <y>  =  [I  +  pR]~*  *  V  (47) 

Equation  (47)  points  out  that  the  control  loops  only  compensate  for  the  average 
of  x(w)  over  the  nulling  band.  This  occurs  because  feedthrough  is  not  sensed 
in  the  covariance  matrix,  R,  as  are,  for  example,  weighted  channel  frequency 
variations  (e.g.,  those  occurring  in  the  front  end  filters).  Using  (45)  and 
(47),  and  noting  <Ay>  =  0,  then  the  adapted  output  power,  Pq,  is  given  by 

po  =  V+  •  [I  +  pR]'1  •  R  •  [I  +  pR]'1  •  V  +  <Ay+  •  R  •  Ay>  (48) 

The  first  term  in  (48)  is  the  output  of  an  ideal  processor  (i.e.,  one 
with  no  feedthrough).  The  second  term  is  a  residual  component  of  power  due 
to  the  frequency  variation  in  the  feedthrough.  This  residual  power  is  weight 
independent,  and  cannot  be  compensated  for  by  the  adaption  process! 

The  above  limitations  are  quite  significant  and  have  considerable  impact 
on  system  performance.  To  illustrate  this  impact,  let  us  rewrite  Eq.  (48) 
in  the  form 

PQ  =  pa (without  feedthrough)  +  k  P^n  (49) 
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That  is,  the  actual  output  of  the  nulling  processor  consists  of  the  theoretical 
adapted  output,  assuming  perfect  weights,  plus  a  residual  term  proportional  to 
the  input  power.  Figure  10  illustrates  the  effect  of  this  latter  term  on  system 
performance.  In  Fig.  10,  we  plot  the  input-output  power  relationship  for  the 
nulling  processor  either  in  the  quiescent  (unadapted)  state  or  the  adapted  state. 
The  difference  between  the  two  curves  is  indicative  of  the  cancellation  realized 
by  the  processor.  Clearly  the  constant  k,  dependent  on  Ay,  limits  the  achiev¬ 
able  cancellation.  In  the  general  case  (i.e.,  non-broadside  incidence),  the 
constant  k  is  dependent  on  the  feedthrough  parameters  (i.e.,  isolation,  time 
delay,  etc.),  interference-source  angle  of  arrival,  and  quiescent  steering 
vector.  For  weight  circuits  designed  to  be  as  identical  as  possible,  k  would 
be  near  maximum  for  broadside  incidence.  Consequently,  it  is  for  this  case 
which  we  estimate  quantitatively  the  cancellation  dependence  on  feedthrough 
in  the  following  paragraphs. 

To  obtain  an  estimate  of  the  constant  k,  consider  the  second  term  in 
Eq.  (48) .  We  have 

'V  feedthrough  '  <*I+  '  5  '  *X>  ‘  »j  £+  •  <tW>  •  H  (50) 

The  cancellation  C  is  defined  in  Eq.  (16).  Assuming  the  feedthrough  term 
dominates  the  adapted  output  power,  then  (I/N)  is  given  by 

R.  E^"  •  <Ay  Ay^  •  E 

d/N)a  -  ~  (51) 

(w  +  <y>)  •  (w  +  <y>) 

In  order  to  simplify  Eq.  (51),  note  that,  for  broadside  incidence,  E^  is  propor¬ 
tional  to  the  first  eigenvector,  ,  of  R.  Furthermore,  assuming  an  earth 
coverage  quiescent  steering  vector,  it  can  be  shown  that  F.q.  (47)  reduces  to 

*  +  <X>  =  colli  -  J  Jj  (52) 

valid  for  broadside  incidence.  Hence 

(w  +  <y>)  +  •  (w  +  <y>)  =  1  -  ^  (55) 
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OUTPUT  POWER,  F^(dB) 


Fig.  10.  Output  power  vs.  input  power  for  a  nulling  processor 
illustrating  the  effect  of  weight  feedthrough  on  the  cancellation 
performance. 
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Consequently,  the  expression  for  (I/N)  becomes 

2L 


(I/N)  =  R. (E  •  E) 

a  J  —  ~ 


+  •  <Ay  Ay  >  •  £1 


1  -  1/N 


t  2 

Now  R.E  •  =  R.  £  |E.  j  =  Trace(R)  =  s.  =  maximum  eigenvalue  of  R.  Hence, 

3  3  k=l  K  =  1 

(I/N)a  =  nTT  si  ®i  *  <£V  £y+>  *  (55] 

To  compute  ( I/N) ^  for  the  same  earth  coverage  quiescent  constraint,  we  note 


(I/N), 


V  •  R  •  V 
V+  •  V 


=  sx/N 


so  the  cancellation  C  is  given  by  (see  Eq.  (16)) 


C  =  Ne^  •  ^A y  Ay+^  * 


Clearly,  the  cancellation  is  directly  dependent  on  the  frequency  variation  of 
the  feedthrough.  Equation  (57)  can  be  estimated  for  two  characteristically 
different  types  of  feedthrough  --  correlated  and  uncorrelated. 

(a)  Uncorrelated  Feedthrough  --  One  example  of  uncorrelated  feed¬ 
through  which  might  occur  is  when  the  feedthrough  path  length  from  weight 
input-output  is  the  same  as  the  weighted  path  length  from  input-output,  but 
the  isolation  ^  Ay  Ay^  takes  the  form  of  a  diagonal  matrix 

r°2  M 

<Ay  Ay+  >  =  |yQ|2  I  o2  (5f 


2  2 
where  o  is  the  isolation  tracking  error,  and  |y  |  the  isolation  level  of  the 

feedthrough  path  relative  to  the  weighted  path.  Then  (57)  reduces  to 

C  %  N|yo| 2  a2  (59) 


42 


which  is  the  result  tabulated  in  Table  1.  Equation  (59)  assumes  an  earth 
coverage  steering  vector.  Since  (I/N)^  is  dependent  on  steering  vector, 
whereas  (I/N)  is  an  absolute  level  independent  of  steering  vector,  Eq.  (57) 

cL 

specifying  C  would  change  dependent  on  V.  Otherwise  said,  referring  to 
Fig.  10,  the  unadapted  curve  in  Fig.  10  moves  up  or  down  dependent  on  V, 
whereas  the  nulled  output  in  the  feedthrough  limited  region  remains  independ¬ 
ent  of  V. 


(b)  Correlated  Feedthrough  —  Correlated  feedthrough  errors  occur 
when  a  time  delay  difference  exists  between  the  weighted  and  feedthrough 
paths.  Effects  on  cancellation  due  to  this  type  of  error  can  be  estimated 
assuming  identical  time  delay  mismatches,  defined  by  t,  for  each  weighting 
circuit.  Then  x.(w)  takes  the  form 

j  IV  (c»i  T) 

Y(w)  =  Xq  e  °  (60) 

It  is  a  straightforward  procedure  to  evaluate <Ay  Ay+^.  We  obtain 


<Ay  Ay+> 


t 


j  AWu)qT 


sin  it  BWt  I  2v 
ttBW  •  T  1  ' 


(61) 


where  BW  =  bandwidth  (Hertz),  and  AW  =  W-l.  For  BW  •  T  <<  1,  we  obtain 

,2 


<Ay  Ayt>  =  ^  ^  L 


BW  »  T) 


(62) 


Substituting  Eq.  (62)  into  (57),  and  noting  that  ej  =  —  [1,1,  ....  1],  we  obtain 

2  1  /N 

CU2  |  Y0 1  2  (lT  -y—  (63) 


Equation  (63)  is  the  result  cited  in  Table  1. 

Finally,  although  not  related  to  the  feedthrough  effect  developed  above, 
it  remains  to  evaluate  the  effects  of  using  quantized  weights,  as  alluded  to 
in  Table  1.  It  is  interesting  to  note  that  analyzing  this  effect  is  similar 
mathematically  to  that  for  the  feedthrough  analysis.  To  see  this,  we  express 
the  error  in  weight  setting  in  the  form 


w  =  w  ♦  Aw 


(64) 
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where  w  is  the  precise  (desired)  weight  and  Aw  is  the  error  due  to  quantization. 
Define  the  rms  quantization  error  according  to 


,2  =  E  j  1 1  -3L  I2j  »  ~2  EiAwk^ 

w.  w, 


a  =  till  — ^  |  >  =  ~2  Awk  ^  (65) 

"k  ’  wk  . 

where  the  expected  value  is  statistically  a  function  of  all  possible  weight 
settings.  The  adapted  interference  to  thermal  noise  ratio  is  given  by 

e{w^  •  R  •  w}  e{aw^  •  R  •  Aw} 


(1/N). 


— t  — 

w  •  w 


E{w+  •  w} 

We  can  express  Eq.  (65)  in  the  same  form  as  (51),  i.e., 

E^  •  e{aw  Aw^  }  •  E 


(66) 


(I/N). 


— t  — 

w  •  w 


(67) 


Assuming  uncorrelated  errors  between  channels,  and  using  Eq.  (65) 

V-2 


s|aw  Aw+ |  = 


'1 


W2 


-2 

WN 


(68) 


Assuming  an  earth  coverage  steering  vector  constraint,  we  have  w^ 
*2  =  •••  ■  WN  =  0-1/N  so  that  (68)  reduces  to 

.2 


«  1-1/N, 


;|aw  Aw+|  = 


(1-1/N) 


1/N 


1/N 


(69) 


Using  Eq.  (69)  in  (67),  we  obtain 


(I/N). 


2  ,,  1>2 
S.O  (1  -  |j) 


l-h 


Sj,a2/N 


(70) 
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Since  (I/N)^  =  s^/N,  it  follows  that 

C  =  a2  (71) 

which  is  the  result  referred  to  in  Table  1. 

Another  class  of  errors  alluded  to  in  Table  1  which  occurs  in  the  weight¬ 
ing  circuits  is  I/Q  imbalance  from  true  quadrature  separation.  Because  of  the 
inherent  compensation  of  these  errors  by  the  action  of  a  feedback  processor, 
discussion  of  this  error  type  is  deferred  to  the  next  section. 

D.  Effects  of  Component  Imperfections  in  the  Feedback  Loops 

The  function  of  an  adaptive  processor  is  to  provide  control  of  the 
weights  following  the  antenna  port  outputs.  Many  techniques  for  adaptive 
control  of  the  weights  have  been  developed  and  used  in  practice.  Each  has 
its  own  peculiar  properties.  Since  the  feedback  type  of  processor  employing 
either  analog  or  digital  processing  is  most  commonly  used,  we  direct  our 
attention  to  it. 

As  discussed  in  Section  III,  one  of  the  basic  functions  of  the  feedback 
loops  is  to  "sense"  the  correlation  matrix  of  the  N  signal  channels.  In 
practice,  the  correlation  matrix  sensed  by  the  control  loops  includes 
imperfections  in  both  the  signal  channels  and  feedback  channels.  If  these 
imperfections  are  frequency  independent,  then  there  exists  a  unique  set  of 
weights  that  minimizes  the  output  power,  therefore  compensating  for  the 
errors.  It  is  not  surprising  that  a  feedback  type  processor  can  adapt  to 
this  unique  set  of  weights  because  the  inherent  nature  of  the  feedback 
mechanism  is  to  drive  the  output  toward  zero.  Consequently,  we  anticipate 
that  the  main  effect  of  frequency  independent  errors  is  their  impact  on  loop 
dynamics.  Otherwise  said,  as  long  as  the  feedback  loops  remain  stable,  the 
feedback  processor  is  considerably  tolerant  of  frequency- independent 
imperfections. 
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To  demonstrate  the  above  conjecture,  and  to  develop  a  quantitative 
assessment  of  component  error  effects  on  loop  dynamics,  we  consider  the  model 
of  the  analog  feedback  processor  illustrated  in  Fig.  11.  The  following  error 
types  will  be  considered: 

1. 

Pre-weighting  channel  tracking  errors:  H^(w) 

2. 

Post-weighting  tracking  errors:  (^(w) 

3. 

Feedback  loop  tracking  errors:  F^Cw) 

4. 

I/Q  mixer  imbalance  errors:  and  £2^ 

5. 

I/Q  weight  imbalance  errors: 

6. 

Loop  gain  imbalance:  ^u/^k 

7. 

Mixer  dc  offsets:  Vj^c  anc*  V2k^C 

As  before,  A^(u)  denotes  the  antenna  frequency  response  for  the  ktl1  channel. 

We  assume  a  single  pole  filter  (F(s)  =  1/(t0s+1))  in  each  feedback  loop.  For 
an  N-channel  system  with  component  errors,  the  equations  characterizing  the 
dynamic  system  response  must  be  formulated  as  a  2N  dimensional  problem.  In 
the  following,  we  first  develop  the  general  set  of  2N  equations  governing  the 
system  response  and  then  consider  separately  the  effects  of  frequency  inde¬ 
pendent  and  frequency  dependent  errors. 

Denote  the  weights  for  the  I  and  Q  channels  as  and  3^,  respectively, 

where  k  =  1,  ...,  N.  In  a  statistical  formulation,  and  3^  are  character¬ 

ized  as  random  processes  which  are  dependent  on  the  interference  noise 
process  incident  on  the  antenna  and  the  time-evolution  of  the  expectations 
E(a^}  and  E{S^}  (denoted  hereafter  as  a^,  3^,  for  simplicity)  is  governed  by 
2N  differential  equations.  Statistical  variations  relative  to  the  expected 
values  of  the  weights  are  called  "weight  noise"  and  are  treated  separately  in 
Section  IV-E.  From  Fig.  11,  we  find  that  the  2N  equations  for  and  3^  are: 

Vk  *  l»k  *  E1Ik>!  ■  vik 

k  ■  1 . n  (72) 

Vk  V  *  »k  *  E(V)  '  V2k 
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S’  *  <frkEkCq*^> 


Fig.  11.  Model  of  a  wideband,  analog  feedback  nulling 
processor  including  the  effects  of  various  component 
imperfections. 


where  the  complex  steering  vector  +  j and  we  assume  ^_(w)  =  0 

(i.e.,  no  weight  feedthrough)  V.  It  can  be  shown  that  (e.g.,  see  the  analysis 
leading  up  to  Eq.  (13)) 


E1Ik)  ■  "ik 

E(Qk)  -  u2k 


+  J£lk 

Re-je  K  (R"  •  w)R 


,,dc 

+  Vlk 


JE 


2k 


Imje  (R"  •  w)k  >  +  p2k  V 


dc 

2k 


(73) 


In  (73),  V^k  and  V2k  are  the  I/Q  channel  beam  steering  voltages  and  R"  is  the 
N  x  N  covariance  matrix  sensed  by  the  correlation  mixers  (Fig.  11),  and  is 
given  by 

B"  =  <Fk  Ek  Cq  E5>  <74> 

k,q  M  H 

where  (to)  = 
weight  w,  defined  by 

3*k 

\  =  “k  +  j3k  e  (75) 

JW'k 

The  factor  e  characterizing  weight  imbalance  errors,  accounts  for 
the  deviations  from  true  quadrature.  Any  absolute  phase  error  in  the  weight 
transfer  function  is  included  in  the  definition  of  Ck(w). 

Equation  (72)  defines  2N  real  equations  specifying  the  time  evolution 
of  and  g^.  This  evolution  is  dominated  by  the  properties  of  R"  defined  in 
(74) .  Observe  that  R"  is  non-hermetian,  indicating  that  the  eigenvalues  of  R" 
will  in  general  be  complex,  leading  to  oscillatory  solutions  for  and  g^. 

Once  the  and  g^  adapt  to  steady  state  (assuming  that  stable  solutions  to 
(72)  exist),  then  the  adapted  output  power,  P  ,  can  be  computed.  Define  the 
N  *  N  signal  channel  covariance  matrix  (including  signal  channel  imperfections) , 
R' ,  to  be 

R'  =  <C,  E,  C*  E*>  (76) 

=kq  '  k  k  q  q'  ’ 


(cu)  Ak(ai).  Also  in  (73),  we  have  introduced  the  equivalent 
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The  output  power  is  then  given  by 

P  =  •  R'  •  w  (77) 

o  —  =  - 

Equation  (77),  when  w  is  obtained  from  Eqs.  (72)  and  (73),  specifies  the 
dependence  of  adapted  output  power  on  loop  component  errors.  To  determine 
w  in  the  steady  state,  substitute  Eq.  (73)  into  (72)  and  set  =  0. 

The  resulting  expressions,  expressed  in  the  vector,  form  are: 


a  +  RejM^  •  R1’  •  w 


=  V  —  V 
-1  -1 


dc 


+  ImjM^  •  .  R"  •  w 


—2  —  —2 


dc 


where  we  have  defined  the  diagonal  matrices 


(78) 


►je 


lk 


‘lk 


k»q 


-JE 


=  P 


2k 


2k 


k,q 


(79) 


6^  ^  is  the  Kronecker  delta  function  and  (V^0)^  =  Plk  Vj^0,  2^^  E  y2k  ^k^' 

Observe  that  the  effect  of  the  mixer  dc  offset  voltages  is  to  modify  the  beam 

steering  voltage.  This  does  not  affect  the  ability  of  the  processor  to  null, 

but  does  affect  the  adapted  radiation  pattern  away  from  the  null.  This  could 

be  a  serious  consideration  for  the  mode  of  operation  of  steering  a  maximum 

gain  beam  to  a  user  close  to  an  interference  source,  but  in  the  earth-coverage 

mode  it  is  not  a  significant  effect.  In  either  case,  digital  calibration 

techniques  can  be  used  to  compensate  for  this  effect.  Consequently,  in  the 

dc  dc 

following  analysis  we  assume  and  to  be  zero. 


s.,  denote 
N 


Closed  form  solutions  to  Eq.  (78)  are  in  general  difficult  to  obtain. 

The  task  can  be  simplified  considerably  by  assuming  that  the  interference 
sources  are  large  relative  to  loop  threshold.  That  is,  if  s^,  . 
the  N  eigenvalues  of  R",  we  assume  that  p  |  s .  |  >>  1  for  each  interference 
source.  In  this  case,  the  a  and  terms  on  the  lefthand  side  of  Eq.  (78)  are 
negligible.  Equation  (78)  can  then  be  rewritten  as  2N  real  equations  for  the 

solution  of  the  vector  X  where 
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X 


(80) 


=  R"  •  w 


Solving  (78)  for  X,  w  can  then  be  determined  by  inverting  (80) . 


w  =  R1 


-1 


(81) 


Substituting  (81)  into  (77)  yields  the  general  dependence  of  output  power  on 
loop  component  errors. 

It  is  convenient  to  classify  the  effects  of  component  errors  into  two 
categories:  "frequency  dependent"  and  "frequency  independent."  Consider 
first  the  latter  category. 

1 .  Frequency  Independent  Component  Error  Effects  on  Cancellation 


On  substituting  the  adapted  solution,  Eq.  (81)  into  (77),  Pq  takes 


the  form 


P  =  X+  •  (R"-1)+  •  R'  •  R"-1 • 
o  —  =  =  = 


(82) 


If  we  assume  all  errors  to  be  frequency  independent,  then  (82)  reduces  to  a 

simpler  form.  For  this  case,  let  Ck(w)  -  e+;,^k  to  represent  the  net  phase 

error  of  the  kth  channel  and  let  F,  (<d)  =  e^k.  Then  (74)  for  R"  can  be 

K  =k  ,  q 

expressed  in  the  form 

F  F 

R"  =  <7=r~  •  C,  E.  E*  C*  >  =  7^  R.’  (83) 

=k,q  ^C^  k  k  q  q  '  Cfc  =k,q 

so  that  R'  and  R"  are  related  by  a  diagonal  matrix,  i.e.. 


R' 


D  •  R' 


(84) 


where  and  =  0  for  i  /  k.  R'  and  R  (where  R  is  defined  at  the 
output  ports  prior  to  weighting)  are  also  related  by  a  diagonal  matrix.  From 
(76), 


j5k  -jCq 

R'  =  e  K  R.  e  9 
k,q  k,q 


(85) 
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where  R,  =  ^E.E*^  ,  which  can  be  written  as 
k ,  q  k  q 


R'  =  C  •  R  •  C 


(86) 


where  =  e  and  =  0  for  k  t  l.  It  can  be  easily  shown  that  the 
eigenvalues  of  R'  are  the  same  as  those  of  R.  Using  Eq.  (86)  in  (82),  we 
obtain 

P  =  [X+  •  (D'V]  •  R''1  •  [D'1  •  X]  (87) 


which  takes  the  same  form  as  the  error-free  solution  having  equivalent  beam 

steering  vector  =  D-^  •  X.  Since  the  R'  has  the  same  eigenvalues  as  R, 

P  can  assume  as  low  a  value  in  this  case  as  for  the  error  free  case.  Conse- 
o 

quently,  we  conclude  that  frequency  independent  errors,  assuming  the  loops 
remain  stable,  do  not  affect  the  nulling  capability  of  the  feedback  processor 
They  do,  however,  affect  the  adapted  radiation  pattern  away  from  the  nulls. 
This  is  the  result  referred  to  in  Table  1.  These  results  have  been  verified 
by  numerous  computer  simulations. 

2.  Frequency  Dependent  Errors 

A  general  analysis  of  the  effects  of  frequency  dependent  loop 
component  errors  is  difficult  to  carry  out.  It  is  clear  that  in  this  case 
R"  and  R'  are  no  longer  simply  related  by  Eq.  (83);  therefore  the  achievable 
cancellation  is  affected  by  frequency  dependent  errors.  Setting  F^tio)  =  1 
and  C^(w)  =  1  leads  to  the  analysis  channel  tracking  errors  (mismatches  in  the 
H^(w))  presented  in  Ref.  1.  Consequently,  the  major  thrust  of  this  Section  is 
to  include  the  effects  of  Fj.(w)  and  C^(w)  on  achievable  cancellation. 

To  accomplish  this,  we  resort  to  the  simple  two-channel  model  illustrated 
in  Fig.  12,  assuming  broadside  incidence  (A^  =  A2) .  For  this  model,  assume 
that  the  interference  source  is  large  relative  to  threshold  (us^  >>  1).  The 
expression  for  the  output  frequency  spectrum,  Vo(w),  can  be  shown  to  be 


/o(w)  “  Vl[<F2C2>Cl(u)  -<F2Cj>  C2(U)J 

+  V2^F1C1>  C2(<i,)  -<FIC2>Cl(u))] 


(88) 
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We  define  the  frequency  response  errors,  AF(w)  and  AC(w),  to  be: 


AF(io)  =  F2(ai)  -  F j (to) 

A C(w)  =  C2(u))  -  C^ (w)  (89) 

Using  Eq.  (88),  we  can  show  that  the  effects  of  AF  on  cancellation  are  of 
secondary  importance  relative  to  AC.  That  is,  differences  between  the 
feedback  loop  frequency  response  functions  have  only  a  weak  effect  on  the 
achievable  cancellation  level  compared  to  differences  in  the  response  func¬ 
tions  of  the  signal  paths.  To  see  this,  substitute  Eq.  (89)  into  (88).  We 
obtain 

V“)  =  V*  [<fJ  AC>  CjOo)  -<F*  Cj>  AC  (to) 

+  <AF*  AC>C1(w)  —  <AF*  Cj>  AC(tu)  ] 

+  V2[<F*  Cj>AC(«)  —  <F*  AC>  Cj  (io)  ]  (90) 


Equation  (90)  characterizes  the  output  frequency  spectrum  of  the  output  signal 
voltages.  For  an  ideal  system,  VQ(w)  would  be  zero.  The  fact  that  it  is 
not  is  attributable  to  the  tracking  errors,  AC(w)  and  AF(u>).  Note  that  all 
terms  involving  AF  are  of  order  AF  •  AC,  indicating  only  a  weak  dependence 
of  cancellation  on  AF  relative  to  its  dependence  on  AC.  Note  that  if  AC  =  0, 
v0(“)  =  0  independent  of  AF.  Consequently,  we  conclude  that  the  dominant 
effect  on  cancellation  is  due  to  those  frequency  dependent  tracking  errors 
in  the  signal  channel  for  which  the  results  of  Ref.  1  apply  directly. 

3.  Effects  on  Loop  Stability 


Consider  now  the  stability  requirements  imposed  by  the  set  of  2N 
differential  equations,  (72)  and  (73).  It  is  convenient  to  recast  these 
equations  into  the  form  of  a  single  2N  x  2N  matrix  equation: 


(91) 
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where 


(Mk,q 


+  j£ 


p  ,  Re<R''  e 
lk  I  =  k , q 


lk 


and 


(^k,q 


(5l>k.q 


^k.q 


Ulk  Re 


j  R"  e 
=  k,q 


j(VClk} 


U-,,  I  m  <  R' ’  e 

2k  J=k,q 


2k 


u2k  Im 


■  jRL* 

-k,< 


(92) 


We  consider  only  the  narrowband  (FBW  «  1)  case  in  which  path  length  errors  in 
the  various  parts  of  the  loops  can  be  approximated  as  fixed  phase  errors. 
Recall  that  R"  implicitly  includes  a  frequency  independent  tracking  error 

'jCq 

C  (a)  =  e  M  due  either  to  weight  circuit  hybrid  mismatches  or  to  signal 
channel  path  length  errors.  If  we  define  the  2N  x  2N  stability  matrix  S  as 

S  = 

where  S  is  a  real  matrix,  then  the  eigenvalues  of  S  determine  the  stability  of 
the  coupled  set  of  2N  loops.  In  particular,  if  we  denote  the  eigenvalues  of 
S  by  X^(£=l  ...,  2N),  then  stable  operation  requires  that 


1%  =2 


=1  i+! 2J 


(95) 


ReU^)  >  0  (94) 

for  {.  =  1,  ...,  2N.  In  general,  S  is  not  a  symmetric  matrix  and  the  A „  are 
in  complex,  admitting  oscillatory  or  unstable  solutions  for  and  Bj.. 

Example;  A  Single  Loop 

A  general  study  of  the  stability  criterion  defined  by  Eq.  (94)  is 
complicated  and  requires  numerical  solutions  for  the  eigenvalues  of  Eq.  (9A) . 
However,  the  analysis  for  a  single  loop  is  analytically  tractable  and  it  is 


instructive  to  examine  this  example  to  obtain  insight  into  the  dependence  of 
stability  on  the  various  component  errors.  A  single  loop  is  shown  in  Fig.  13. 
If  it  is  stable,  it  will  simply  tend  to  turn  weight  "off"  when  a  signal  is 
applied  to  the  input. 

For  the  case  of  a  single  loop,  we  adopt  the  following  notations 


(1) 

R 

0 

input  power  to  loop  relative  to  thermal  noise 

(2) 

U11  = 

v,  =  nn,  where  n  =  gain  mismatch  between 

I  and  Q  channels 

(3) 

£11  = 

E 1  *  £2i  =  £2*  ^l  =  ^ 

(4) 

C(uj)  = 

-j(A_-A ,)  .. 

-]5  r,  ,  J  2  1  .  -jA 

e  J  F(oi)  =  e  =  e 

With  these  definitions,  the  stability  matrix  S  takes  the  form 


S  = 


1  +  y  R  cos  e ' , 
o  1 


nu  Rq  sin  e'2 


-  u  RQ  sinO+e  ' 


1  +  nn  Rq  cos  ('I'+e  '  2) 


(95) 


where  e'j  =  e j  -  (A-£),  e'2  =  £2  *  (A-£) .  The  two  eigenvalues  of  S  can  be 
computed  directly,  and  are: 

A 1  2  =  [1  +  j  u  Rj(cos  e’j  +  n  cos(i^  +  e'2))] 

12  1 

j  cos  e'j  -  j  n  cos  e'^  cos(i|>  +  e'2) 

+  j  cos^(i);  +  e'2)  -  n  sine'2  sinO  +  e’^)  ^ ^  (96) 

Applying  the  stability  criterion  of  Eq.  (94),  Eq.  (96)  leads  to  the  stability 
diagram  illustrated  in  Fig.  14.  In  this  diagram,  the  errors  n,  <1*  and  (£-A)  are 
treated  as  parameters,  and  stability  is  plotted  as  a  function  of  correlator 


t 


A=A|-A2 


€2*s€2  +  ( 


Fig.  13.  Loop  parameters  used  for  the  single  loop  stability 
analysis  leading  to  the  stability  diagram  of  Fig.  14. 
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Fig.  14.  Stability  diagram  for  a  single,  first-order, 
Applebaum- Howe 11s  type  control  loop. 


output  phase  errors  e^  and  The  loop  is  stable  if  e^,  £2  lie  inside  the 

solid  boundary  indicated  in  the  Figure.  The  boundary  is  determined  from  the 
two  conditions. 

cos  c'j  +  n  cos(t|;  +  e^)  >  0 

W2  -  e’jl  <  90° 

which  follow  directly  from  (94)  and  (96),  assuming  y  Rq  >>  1 
(98)  we  have  defined 


e'j  =  Ej  +  U-A) 

(99) 

£ '  2  =  e2  + 

(100) 

indicating  that  differential  phase  errors  in  the  feedback  paths  are  equivalent 
to  identical  correlator  phase  offsets  in  the  I  and  Q  channels.  Generally 
speaking  there  is  a  biased  region  of  stability  about  eJ  =  =  0, 

To  obtain  a  physical  feel  for  how  the  loops  become  unstable,  consider 
a  simple  unstable  condition  and  track  through  the  feedback  process.  Basically, 
stability  (or  instability)  is  a  consequence  of  the  negative  (or  positive) 
feedback  condition  existing  at  the  output  of  the  correlator,  where  errors  in 
the  mixer  generate  a  drive  voltage  in  the  quadrature  channel.  Once  the 
quadrature  channel  is  turned  on,  the  output  now  has  a  component  shifted  90°. 
Hence,  should  the  correlator  differential  error  between  I/Q  channels  exceed 
90°,  this  90°  shift  combines  with  the  quadrature  weight  90°  shift  resulting 
in  a  180°  phase  reversal  at  the  sum  junction  --  i.e.,  the  feedback  in  the 
quadrature  channel  is  now  positive  feedback,  resulting  in  an  increasingly 
large  growth  in  signal  level  driving  the  weight,  and  the  loop  becomes 
unstable.  Thus,  the  stability  boundary  as  a  function  of  e'j,  e'2  inter¬ 
sects  ±  90°  on  all  axes  of  the  diagram.  When  £j  =  -e7,  a  45°  error  in  both 
1  and  Q  yields  a  90°  differential  output  error  in  the  correlator  and  an 
unstable  condition,  as  can  be  seen  on  the  stability  diagram.  Errors  in  the 
hybrid  (ip)  or  differential  path  length  (t-A)  result  in  a  more  complicated 


(97) 

(98) 

In  (97)  and 
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I/Q  coupling,  but  essentially  lead  to  some  type  of  positive  feedback  when  the 
loops  become  unstable. 

Some  N-Loop  Examples 

The  above  results  indicate  that  the  stability  characteristics  of  a 
single  loop  are  quite  tolerant  to  most  reasonable  values  of  component  errors. 
Practically  speaking,  one  would  anticipate  that  and  e'2  would  be  5  to  10° 
at  most  and  that  weight  hybrid  errors,  would  be  less  than  1  to  2°.  Loop 
gains  might  differ  by  as  much  as  2  or  3  dB.  Consequently,  it  appears  to  be 
relatively  easy  to  design  a  single  stable  loop.  The  question  then  arises, 

"How  does  the  processor  stability  depend  on  the  number  of  loops  and  number 
of  interference  sources?" 

It  is  useful  to  consider  a  simple  special  case  in  order  to  obtain  a  feel 
for  the  dependence  on  the  number  of  loops.  Consider  the  N-loop  problem  having 
equal  I  and  Q  mixer  imbalance  errors  in  each  loop.  Denote  this  error  by  e. 

Then  we  have 

t  a  +  a  +  y  Re[eJE  R'  •  wl  =  V, 

o—  —  =  —  1 

toA  +  !  +  P  Im[e^e  R'  •  w]  =  V2  (101) 

Equation  (101)  can  be  rewritten  in  complex  notation  as 

t^w  +  w  +  neJe  R'  •  w  =  Vj  +  jV2  (102) 

where  now  w  =  a  +  j&  since  we  assume  ^  =  0.  Denote  the  eigenvalues  of 
R'  by  s^.  Then  the  stability  of  (102)  is  governed  by  the  eigenvalues 
=  1  +  us^  eJE.  Clearly,  if  Re(A^)  -  0,  we  must  require  that  e  <  90°. 

This  is  the  same  stability  requirement  as  for  a  single  loop.  Consequently, 
we  anticipate  that  increasing  the  number  of  loops  should  not  adversely  affect 
the  stability  of  the  processor. 
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In  practice,  the  various  component  errors  will  be  statistically  distributed, 
with  some  channels  having  little  error  and  others,  perhaps,  a  larger  error. 

The  above  result  then  is  a  worst-case  condition  for  this  more  realistic  case 
and  we  anticipate  the  stability  of  the  processor  will  improve  as  the  number 
of  loops  increases.  To  examine  this,  we  resort  to  a  numerical  simulation  of 
Eq.  (93)  including  all  types  of  component  errors.  Each  component  error  will 
be  assumed  to  be  uniformly  distributed  over  a  range  within  a  maximum  error 
bound.  Each  simulation  trial  uses  a  set  of  errors  randomly  selected  from  each 
distribution.  The  stability  criterion  for  each  sample  set  of  errors  can  be 
checked  by  computing  the  eigenvalues  of  S  and  applying  the  condition 
Re(A?)  ■>  0,  H=1 ,  ...,  2N.  This  process  is  repeated  for  a  number  of  trials 
with  a  specified  error  bound  for  each  component.  The  percentage  of  cases 
for  which  the  processor  is  unstable  is  tabulated.  In  this  manner,  the  proba¬ 
bility  of  achieving  stability  can  be  studied  as  a  function  of  the  error 
bounds  for  a  particular  component.  The  result  can  then  be  used  in  designing 
the  processor  to  a  fixed  set  of  specifications. 

This  procedure  can  be  utilized  to  examine  processor  stability  as  a 
function  of  antenna  type,  number  of  loops  and  number  of  interference  sources. 
Consider  the  results  of  Fig.  15,  where  we  show  processor  stability  as  a 
function  of  the  number  of  loops  for  a  single  interference  source.  For 
example  purposes,  we  use  the  7-element  hexagonal  array  antenna  geometry 
illustrated  in  the  figure,  and  examine  the  percent  of  unstable  eigenvalues 
vs  error  bound  as  the  number  of  control  loops  increases  from  1  to  7.  For 
simplicity,  we  consider  only  I/Q  mixer  imbalance  errors  to  illustrate  the 
dependence  ou  N.  Observe  that  for  a  single  loop,  no  value  of  mixer  imbalance 
errors  leads  to  an  unstable  loop  until  the  e  exceeds  45°.  This  critical 
value  allows  the  case  Ej  ^  -  e,  to  occur  for  some  set  of  random  errors.  As 
mentioned  previously,  this  results  in  an  unstable  loop.  More  important ly , 
observe  in  Fig.  15  that  as  the  number  of  loops  increases,  the  loop  stability 
becomes  tolerant  of  a  larger  error  bound,  as  predicted  above.  Clearly,  for 
7  loops,  the  tolerances  on  I/Q  mixer  imbalance  errors  would  be  very  easy  to 
sat i sfy . 
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Fig.  15.  Stability  dependence  as  a  function  of  I/Q  mixer 
imbalance  errors  for  1,  2,  4  and  7  control  loops.  Antenna 
geometry  used  was  7-element  hexagonal  array  haveing  D/X  =6.6 
with  a  single  interference  source  positioned  at  9  =  4°,  <f>  =  0° 
FBW  =  0. 
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We  now  fix  the  number  of  loops,  and  vary  the  number  of  interference 
sources.  Consider  the  same  7-element  hexagonal  array  with  6  sources  located 
on  the  FOV  as  shown  in  Table  3  below. 


TABLE  3 

RANDOM  LOCATIONS  OF  SIX  INTERFERENCE  SOURCES 


1  1°  -  10° 

2  4°  30° 

3  8°  80° 

4  5°  135° 

5  7°  205° 

6  6°  280° 


The  probability  of  achieving  stability  as  a  function  of  the  error  bound  on 
mixer  imbalance  errors  is  indicated  in  Fig.  16  for  J=l,  4  and  6.  Clearly, 
the  sensitivity  to  component  tolerance  errors  is  most  severe  for  the  case  of 
six  sources  incident  on  the  7-element  array.  In  this  case  the  processor  has 
used  all  its  degrees-of-freedom  to  null  sources.  Note,  however,  as  suggested 
by  our  simple  example  above,  the  allowable  component  tolerances  never  become 
more  severe  than  those  for  a  single  loop. 

Finally,  consider  a  7-element  array  with  a  double  triangle  configuration 
(Fig.  17).  Again,  we  use  the  6  interference  sources  of  Table  3.  All  the 
errors  (Cj^,  and  A^)  are  allowed  to  vary  statistically  as  in  the 

previous  simulations.  In  this  case,  we  tabulate  the  statistics  of  the  trial 
cases  in  the  following  manner:  We  treat  n  (gain  mismatch)  as  a  parameter 

and  generate  a  set  and  A  k=l . 7.  Then  for  this  random  set  of 

component  errors,  we  generate  a  subset  of  and  c^j,,  k=l,  ...,  7  bounded  by 
c  ,  .  We  consider  twenty  sample  cases  for  e  .  and  c,.  for  each  set  vl<,  ,  \ . 
and  A^.  We  then  store  the  percent  of  unstable  values  as  a  function  of 
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Fig.  16.  Stability  dependence  as  a  function  of  I/Q  mixer 
imbalance  errors  for  1,  4  or  6  interference  sources  of 
Table  3  incident  on  a  7-element  hexagonal  array.  D/X  =  6.6, 
FBW  =  0. 
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<MAX  *  MAXIMUM  i/q  mixer  imbalance  error 

Fig.  17.  Stability  dependence  vs.  I/Q  mixer  imbalance  errors  when  all 

component  parameters  e.  ,  e7  <1/  ,  £,  and  A,,  k  =  1 . 7  vary 

^v-  i-  ^  K 


k’  k 

statistically  using  _ ,  E  „  and  A _  =  20°.  Antenna  configuration 

'  b  max  max  max  y  • 

used  was  7-element  double  triangle  array  having  D/X  =  18.4.  The  six 

interference  sources  of  Table  1  were  assumed  incident  on  the  array. 

FBW  =  0.  Results  are  for  n  =  1  and  n  =  .7. 
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e  ,  and  repeat  for  a  new  set  of  ,,  ,  '  and  .  for  each  set  of  ,  E, 
max  r  k  K  K  k  k 

and  A^,  we  store  the  cumulative  maximum  of  the  percent  unstable  values  obtained 

vs  e  .  Figure  17  illustrates  this  cumulative  maximum  (i.e.,  the  worst  case 
max 

result)  vs  e  using  ip  =  5  =  A  =  20°  for  n  =  1  and  n  =  0.7.  The 

1  max  max  max  max 

results  are  essentially  the  same  for  n  =  1.0  and  n  =  0.7.  Observe  that 

the  curve  of  Fig.  17  has  shifted  only  slightly  to  the  left  of  that  for  Fig.  16, 

"or  which  case  £  =  A  =  tp  =0.  Consequently,  as  predicted  by  Fig.  14 

nia.x  max  max 

■  r  a  single  loop,  errors  in  A^  and  are  not  cumulative  relative  to 
and  e^-  *n  fact>  as  Fig-  14  and  the  above  simulations  indicate,  the  loops 
are  more  tolerant  to  errors  in  ^  and  A^  than  for  and  c^.  We  conclude 
therefore  that  even  loose  tolerances  on  error  parameters  still  allow  stable  loop 
operation,  even  in  the  worse  case  of  N-l  sources  incident  on  an  N-element  array. 

E.  Weight  Noise 

Heretofore  when  the  time  variation  of  the  weights  has  been  mentioned  we 
were  speaking  of  the  statistical  or  ensemble  average  of  the  weights.  It  was 
pointed  out  in  Section  II,  Background,  that  the  differential  equations  for 
w(t)  were  derived  in  terms  of  a  statistical  average.  In  this  section  we  consider 
the  statistical  fluctuations  of  w(t)  around  its  average  value.  We  refer  to 
these  fluctuations  as  weight  noise.  In  an  analog  processor  weight  noise  may 
arise  from  two  sources.  The  first  is  thermal  noise  in  the  weight  control 
circuits.  Secondly,  if  the  input  signals  are  noisy,  some  of  this  noise  is  fed 
back  to  the  weights,  causing  them  to  fluctuate.  In  digital  systems,  quantiza¬ 
tion  and  round-off  errors  can  contribute  to  weight  noise. 

Weight  noise  has  two  observable  effects.  First,  it  increases  the  adapted 
output  power  of  the  processor  over  what  it  would  be  if  the  weights  were  fixed 
at  their  average  value.  Second,  even  when  the  nulling  processor  is  in  a  steady- 
state,  adapted  condition,  the  continuous  small  weight  fluctuations  make  it 
appear  to  be  a  time-varying,  multi-path  channel  to  the  communications  signals. 
Thus  it  will  tend  to  smear  or  broaden  the  signal  spectra.  This  reduces  the 
detectability  of  a  signal  and  can  also  cause  cross-talk  between  signals  that 
are  separated  in  frequency.  In  this  section  we  will  quantify  these  two  effects. 
We  will  calculate  the  weight  noise  contribution  to  total  output  power  and  find 
that  it  is  insignificant  as  long  as  the  processor  dynamic  range  is  not  exceeded. 
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1.  Circuit  Noise 


Circuit  noise  effects  can  be  modeled  as  in  Fig.  18(a).  We  will 

assume  the  input  v(t)  to  be  a  noiseless  process  (i.e.,  a  sinusoid)  so  that 

circuit  noise  is  the  only  source  of  weight  fluctuation.  The  total  weight 

vector  is  the  sum  of  a  noiseless  average  value  w(t)  and  a  noise  vector  n(t) 

whose  components  are  identical,  but  statistically  independent,  zero  mean  random 

~ 

processes  of  variance,  n  .  An  example  of  such  a  noise  source  in  a  baseband 
feedback  loop  is  the  1/f  noise  in  the  video  amplifier  following  the  correlation 
mixer. 

Total  Output  Power: 

Assume  the  processor  to  be  in  a  steady-state  adapted  condition 
for  which  w(t)  is  a  constant.  Consider  the  total  output  power 

Po  =  E(vo(t)  v*(t)}  (103) 

where 

vQ(t)  =  (w+  +  n+(t))  •  v(t)  (104) 

By  substituting  (104)  into  (103),  and  by  using  the  statistical  independence 
of  the  noise  generators,  it  can  easily  be  shown  that 

p  =  w+  •  R  •  w  +  n2  P.  (105) 

o  —  =  —  m  ’ 

•J* 

where  R  is  the  correlation  matrix  E{v(t)v(t)  },  w  •  R  ♦  w  is  the  noiseless 
adapted  output  power  and  P^n  is  the  total  input  power  to  all  channels. 

When  testing  a  processor  one  often  measures  the  dependence  of  PQut  on 
P^n  for  a  single  interference  source,  usually  at  broadside  incidence.  When  a 
narrow-band  test  signal  is  used,  weight  noise  manifests  itself  as  shown  in 
Fig.  18(b)  (for  a  hard-limited  processor).  When  P^n  exceeds  a  particular  value 
(40  dB  in  this  case),  the  weight  noise  dominates  the  output  power  and  Pq 
becomes  proportional  to  P^n  as  predicted  by  Eq.  (105).  Measurements  to  this 
effect  will  be  presented  in  Section  V. 


RELATIVE  OUTPUT  POWER  Idb) 


RELATIVE  INPUT  POWER  (d8)  FREQUENCY  (MH z) 

lb)  (C) 

WEIGHT  NOISE 

Fig.  18.  Effect  of  noise  in  the  feedback  loop  circuitry  on  nulling 
system  performance  for  a  hard- limited  processor,  (a)  Circuit  noise 
is  modelled  by  identical  but  statistically  independent  noise  gen¬ 
erators.  (b)  The  effect  of  weight  noise  on  the  adapted  Pour  vs.  Pjn 
characteristic  of  a  processor  with  single,  broadside  incident  source, 
(c)  Effect  of  weight  noise  on  the  output  spectrum  for  a  sinusoidal 
input  of  varying  power  levels,  (d)  The  equivalent  linear  system  for 
weight  noise  arising  from  the  interference  signal. 
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Output  Spectrum: 

Weight  noise  modulates  any  signal  passing  through  the  processor 
and  broadens  its  spectrum.  The  power  spectral  density  of  the  output  is 

SQ(<o)  =  l"  Ro(x)e'jwT  di  (106) 

_  CO 

where  Rq(t)  is  the  output  autocorrelation  function, 

Rq(w)  =  Efvo(t)  v* ( t- x)  )  (107) 

As  in  Section  II,  we  define  a  correlation  matrix  for  the  channel  signal  vector 
v(t)  as  follows: 

[Rv(t)J  =  E(v(t)  v+(t-x)}  (108) 

The  cross-power  spectral  density  matrix  is 

[§„(«)]  =  /  Rv(t)]  e'ja)T  dr  (109) 

With  these  definitions  it  follows  from  Eq.  (104)  that  the  output  power  spectral 
density  is 

S  (to)  =  w+  •  S  (<o)  •  w  +  Tr(S*(u))©S  (w) }  (110) 

o  —  =v  —  n  =v 

where  Tr  is  the  matrix  trace  (sum  of  the  diagonal  elements),  Sn(io)  is  the  weight 
noise  spectrum  and  ©  denotes  the  convolution  integral  divided  by  2tt.  For  the 
case  where  the  signal  vector  v(t)  is  a  tone  of  amplitude  V,  Pin  =  N | V |  and  we 
have 

S  (to)  =  2  it  { I V 1 2  w+  •  w  6(w)  +  P.  S  (u>) }  (111) 

The  first  term  is  the  line  spectrum  of  the  signal  at  the  output  (6 (to)  is  an 
impulse  in  frequency) .  The  second  is  the  weight  noise  contribution. 

As  an  example  of  Eq.  (Ill),  Fig.  18(c)  shows  the  measured  output  spectrum 
for  a  hard-limited  processor.  As  P^n  is  increased,  the  signal  component  is 
clamped  at  a  fixed  amplitude  by  the  adaptive  processor,  but  the  weight  noise 
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spectrum  increases  in  proportion  to  .  Of  course  a  pure  tone  is  not  a 
realistic  interference  signal;  however,  this  test  is  useful  in  revealing 
directly  the  circuit  noise  spectrum. 

2 .  Interference  Noise 

We  will  now  consider  the  weight  noise  generated  by  a  noisy  input,  v(t). 
We  will  first  analyze  a  standard  Applebaum-Howells  loop,  i  e. ,  one  without  the 
hard  limiter  of  Fig.  18(a).  We  will  calculate  its  effect  on  a  tone  of  constant 
amplitude,  V.,  whose  power  is  insignificant  compared  to  the  power  in  v(t).  We 
will  also  find  the  weight  noise  contribution  to  total  output  power  and  show  that 
it  is  insignificant  as  long  as  the  processor  dynamic  range  is  not  exceeded. 
Section  3  extends  this  analysis  to  hard  limited  Applebaum-Howel 1  loops  under  the 
assumption  of  broadside  incident  interference. 

Weight  Noise  Spectrum: 

Let  the  total  weight  vector  be  the  sum  of  its  expected  value, 
w(t),  and  a  zero-mean  weight  noise  vector  Z(t) .  Referring  to  Section  II, 

Eq.  (3),  the  equation  governing  the  total  weight  vector  is 

tq(w  +  Z^)  +  [I  +  p^v  v+]  •  (w  ■*•  £)  =  V  (112) 

where  we  understand  that  w,  Z  and  v  are  functions  of  time.  Let  v(t)  be  a  zero- 
mean,  gaussian,  band-limited  process  of  bandwidth,  BW.  The  basic  assumption 
underlying  the  analysis  of  nulling  processors  is  that  and  u  are  chosen  so 
that  the  closed-loop  bandwidth,  BW^,,  of  the  processor  is  a  small  percentage 
of  the  nulling  bandwidth,  BW.  Consequently,  Z(t)  is  quasistatic  in  comparison 
with  v(t)  and  can  be  treated  as  statistically  independent  of  the  instantaneous 
value  of  v(t)  (although  it  is  not  independent  of  the  past  history  of  v(t)). 

With  this  assumption,  the  expectation  of  Eq.  (112)  yields  the  well-known 
equation  for  w(t) , 

Tq w  +  [I  +  m R ]  •  w  =  V  (113) 

where  R  i  |v(o)  in  Eq.  (108).  Subtracting  Eq.  (113)  from  (112)  yields  an 
equation  for  Z^t) . 

t ^Z  +  [I  +  uR]  •  Z  +  p[v  v'  -  Rj  ■  Z  =  —  m[v  v+  —  R]  •  w  (114) 
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This  is  a  linear,  vector  differential  equation  with  a  time-varying  coefficient. 
Figure  18(d)  is  the  linear  system  representation  of  Eq.  (114).  The  forcing 
function  is 


£(t)  =  -  p[v  v+  -  R]  •  w 


(115) 


In  Fig.  18(d)  we  see  that  the  time-varying  coefficient  has  a  negligible  effect 
compared  with  g(t)  if  |zj  <<  |w|.  We  are  therefore  justified  in  neglecting 
it  under  normal  operating  conditions.  The  differential  equation  for  weight 
noise  is  then 

Mg]  •  z  =  g  (116) 

Now  we  can  obtain  the  power  spectral  density  of  2.  In  multi-dimensional  linear 
systems  the  relationship  between  power  spectral  densities  at  the  inputs  and 
outputs  is 

S  fu)  =  H+(w)  S.  (w)  H(uj)  (11?) 

=out  ='=in= 

where  H(m)  is  the  matrix  of  system  input-output  transfer  functions.  From 
Eq.  (116)  we  have 

S_(w)  =  [(1  +  jwx  )  I  +  wR]  1  •  S  (<o)  •  [(1  -  jwr  )  I  +  pR]’1  (118) 

L.  O  —  -  -g  U  - 

Now,  since  BW  >>  BW^.,  §g(“)  is  a  very  wide  spectrum  compared  to  the  system 
bandwidth  and  we  can  approximate  it  by  a  constant  spectral  density, 

Sg(u))  =  WE[&&+] 

From  Eq.  (115)  we  find  an  expression  for  the  expectation, 

*  2  *  *  *  2  t 

E(gk  g^J  =  u  [  w.  w  c(vk  vi  v  v  }  -  v  (R  •  ww  •  R}  u  (119) 

ij  J 

where  we  have  used  the  fact  that  E{v  v+)  =  R.  The  fourth  moment  in  (119)  can 
be  simplified  using  the  following  identity  for  gaussian,  complex  random 
processes  (Ref.  5,  Eq.  33) 

EfXjX^XjX*}  =  EUjX*}  •  E{X3X*}  +  E { X j X * }  •  E{X*X3)  (120) 
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Using  this  in  Eq.  (119)  yields 


E{£  £+}  =  p2(w+  •  R  •  w)  R  (121) 

(w+  •  R  •  w)  is  the  adapted  output  power,  P  ,  in  the  absence  of  weight  noise. 
From  (118)  we  get 

p2(w+Rw) 

§z(w)  =  - —  [(1  +  j‘i)T0)  I  +  pBI  g [ C 1  -  iUJT0)  I  + 

(122) 

which  is  the  desired  expression  for  the  weight  noise  cross  power  spectral 
density.  Since  the  off-diagonal  terms  of  Sz(id)  are  non-zero,  the  weight 
noise  is  correlated  from  channel  to  channel  unlike  the  case  of  circuit  noise. 

Weight  Noise  Interaction  with  the  Signal: 

Now  let  the  input  contain  a  communications  signal  waveform, 

V  (t),  buried  in  the  interference,  v(t) .  V  (t)  is  sufficiently  small  so  as 
not  to  affect  weight  adaption  or  weight  noise.  Let  us  find  the  power  spectral 
density  of  the  processor  output,  vQ(t).  By  the  same  process  that  yielded 
Eq.  (110),  we  can  show  that  it  is, 

S  (w)  =  w+  •  S  (w)  •  w  +  y  [S„  (w)  ]  .  .  ©  [S  (w)  ]  .  .  +  w^  S  (w)  w  + 

o  —  =s  —  .  t  -Z  J  n  l  =  s  li  —  =v  — 


.1  [S_(u)] .  .  ©  [S  (w)] . .  (125) 

lj  l=Zv  L=v  Jij 

where  §s(u))  is  the  cross-spectral  density  of  V  (t)  and  we  have  assumed  Vg(t) 
and  v(t)  to  be  independent.  The  first  and  third  terms  are  the  signal  and 
interference  throughput  respectively  and  the  second  and  fourth  terms  represent 
weight  noise  modulation  of  the  signal  and  the  interference  respectively.  Let 
us  examine  the  signal  modulation  effect  for  a  sinusoidal  signal.  Under  this 
assumption  the  third  term  in  Eq.  (125)  reduces  to 

s'(co)  =  (^  •  Sz(co)  •  Vs)  (124) 
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To  a  first  approximation,  SjCw)  affects  the  communications  signals  like  additive 
gaussian  noise.  An  important  question  is  how  much  of  the  power  of  Sj(w)  falls 
in  the  signal  bandwidth.  In  frequency  division  multiplexed  systems,  a  more 
important  question  is  how  much  power  falls  into  adjacent  frequency  slots  causing 
cross  talk  between  channels.  Given  the  parameters  of  the  communications  system, 
these  and  other  questions  can  be  answered  once  SJ(w)  is  known. 
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Increase  in  Total  Output  Power  Due  to  Weight  Noise: 


Let  us  return  to  Eq.  (123)  and  consider  the  last  two  terms  which 
are  the  interference  noise  throughput  and  an  additional  component  due  to 
weight  noise  modulating  the  interference.  We  will  show  that  the  latter  effect 
causes  a  negligible  increase  in  power  as  long  as  the  input  power  is  within  the 
dynamic  range  of  the  processor.  In  the  time  domain  the  terms  of  (123) 
characterizing  the  interference  output  power  can  be  written  as 

Po  =  w+-  R  •  w  +  E{Z+  •  R  ♦  Z)  (135) 


where  we  assume  1_  and  w  are  uncorrelated.  Applying  the  linear  transformation 
of  (125)  we  get 

+  N 

e(Z  •  R  •  Z)  =  l  E{  | Z .  |  }  s.  (136) 

i=l  1  1 


i  a  2 

We  obtain  E {  Z .  }  by  integrating  Eq. 

1 

this,  Eq.  (136)  becomes 


P 

o 


,[l  + 


■V 


2t  BW 
o 


(127)  over  u  and  dividing  by  2  it . 


N  s2  -I 

1  =  1  T  +  yS  J 


With 


(137) 


The  result  in  Eq.  (137)  is  similar  to  that  developed  by  Brennan  et  al.^  for 
sampled  data  systems.  The  bracketed  term  is  the  contribution  of  weight  noise 
to  the  total  output  power.  For  a  single  interference  source  we  have 


P 

o 


r. 

a  -  Pin/pT)2  L 


1 

+  2x  BW 
o 


(Pin/PT)2  1 

(1  *  Pin/PT)J 


(138)) 


For  the  weight  noise  to  be  significant  the  input  power  must  satisfy 
Pin/PT  1  2toBW 


(1391 


However,  this  violates  our  basic  assumption  that  the  weights  are  statistically 
independent  of  the  input,  v(t) .  This  assumption  was  necessary  to  get  Eq.  (113). 


It  requires  that  the  weights  remain  essentially  constant  over  an  interval  of 
several  times  the  correlation  time  of  v(t) ,  which  in  turn  implies 


MS 


max 


(—) 

W 


(140) 


because  t  /us  is  the  minimum  time  constant  in  the  processor  transfer  func- 
o  max  r 

tion  (Eq.  122).  Using  s  %  P .  .  we  have  as  the  condition  for  validy  of 
v  M  6  max  in  1 

this  analysis 


P. 

^  <<  t  •  BW  (141) 

P^  o 


which  contradicts  (139).  (In  fact,  a  common  definition  of  the  processor 
dynamic  range  is  P^  /P  <_  (1/5)  BW  *  tq.)  We  conclude,  therefore,  that  we  can 
never  observe  the  effect  of  Z^t)  on  the  total  power  output  if  the  processor 
is  operating  within  its  dynamic  range. 

3 .  Hard-Limited  Processors 

This  analysis  can  be  extended  only  in  a  limited  way  to  the  hard- 
limited  Applebaum -Howells  loop  (Fig.  18(a)).  This  type  of  processor  is 
described  by  Gabriel  (Ref.  2).  The  reason  for  hard  limiting  is  to  decrease 
the  dynamic  range  requirements  on  the  correlation  mixer  and  subsequent  circuits. 
As  mentioned  in  Section  V  of  this  report  and  in  Ref.  7,  hard  limiters  do  have 
a  harmful  effect  on  system  performance  in  that  they  make  it  difficult  for  a 
processor  to  detect  and  null  the  smaller  of  two  signals  of  widely  different 
power  levels.  Because  hard  limiting  is  used  so  often  in  practice,  weight  noise 
effects  in  hard-limited  processors  is  of  interest.  Therefore  the  necessary 
modifications  to  the  preceding  analysis  are  pointed  out. 


In  a  hard-limited  system  the  correlation  mixer  output  is  described  by 


k(t)  v*(t) 
k(t)  ■  -[-vk(t)-l 


(142) 
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I 
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Consequently  the  basic  differential  equation  is 


To(w  +  Z)  +  [I  +  phM](w  +  Z)  =  V 


where  M  is  given  by 


V  .  V  . 

[M]  .  .  =  t^4 

=  ij  I v. ! 


This  can  be  compared  with  Eq.  (112)  for  the  standard  system.  Brennan  and  Reed 
calculated  the  mean  value  of  M  for  a  gaussian  input  (Ref.  6)  and  found  it  to  be 


=  V T 


where  v  is  the  input  power  per  channel  and  R  is  the  channel  signal  correlation 
matrix  E{v  v+}.  (This  result  is  valid  only  for  array  antennas  with  equal  gain 
elements  because  the  power  is  assumed  to  be  equal  in  all  channels.)  With  this 
result,  the  analysis  proceeds  as  before.  The  equation  for  weight  noise  is 


I  +  \l  +  g]  •  Z  =  £ 

L  J 


L  V^2  "j 

The  equivalent  expression  to  Eq.  (119)  is 


t^{£  £  }] k£  =  I.  ^w  E  “j— 

J  {  1  k 


vbv-  v.v„ 

k  i  l  l 


-  y2h2[M  •  W  W+  •  M] 


The  expectation  cannot  be  simplfied  in  general.  However,  if  we  restrict  v(tl 

to  broadside  incidence  (i.e.,  all  channel  signals  are  equal),  the  expecta- 

2 

tion  reduces  to  E(|v|  }.  Equation  (146)  reduces  to 

2  2 

E(g  g+)  =  (1  -  J)  (w+  •  R  •  w)  R  (147) 

v 

which  is  very  similar  to  Eq.  (121)  for  the  standard  system  except  that  it  is 
restricted  to  the  case  of  broadside  incidence.  Proceeding  as  before  the  weight 
noise  spectrum  (equivalent  to  Eq.  (127))  becomes 
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where 


[Sz(w)]n 


y 2h2N ( 1  -  tt/4)  _ 1 _ 

BW  r -  -  2  ,  2  2  ^ 

(1  +  uh/NSj)  (1  +  w  t ^ ) 


t  ^  =  tq/(1  +  uh/Ns^ ) 


(For  broadside  incidence  the  only  significant  eigenvalue  is  s^  =  . )  The 

expression  for  (w^  •  R  •  w)  in  a  hard-limited  system  under  these  conditions  is 


(w  •  R  •  w) 


IV  Pln 
(1  .  /p-/pTr 


which  is  comparable  to  (131).  (Ve  have  defined  the  threshold  for  a  hard- 
limited  system  to  be 


2  2  - 1 

PT  =  (u  h  N) 

The  weight  noise  modulation  of  a  tone  signal  yields  the  output  spectrum 


So^ 


a  I  IV2  (Pin/PT)  x 


BWfl  +  /P~7pI) 

in  T 


f  i  2  2. 

(1  +  U  Tj) 


This  can  be  compared  to  (132)  for  the  standard  system.  The  maximum  value  of 

the  peak  of  the  spectrum  again  occurs  at  P.  =  P^  and  is 

m  T 


(1  -  m/4)  | V  | 

S0(0)  -  - T6~BW -  volts  /hz 


which  is  very  near  the  result  for  the  standard  system.  The  total  spectrum  power 

reaches  a  maximum  at  P.  =4  P^  and  is 

in  T 


~  /  S'(lu)  dw  < 
2tt  ;  ov  J  — 


2(1  -  m/4)  |Vs 

27  r  r_BF 
o 


volts 


which  is  about  6  dB  less  than  for  the  standard  system.  This  is  because  of  the 
narrower  bandwidth  of  a  hard-limited  system  for  the  same  filter  time  constant, 


1 


The  equivalent  to  Eq.  (139)  for  total  output  power  is 


iy  nyy 

(i  .  Vp— 7p7>2 


i  + 


(1  -  */4) 

2t  •  BW 
o 


'Pi„/PT> 


a  *  ^y). 


(154) 


The  weight  noise  contribution  is  therefore  insignificant  until  exceeds 
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which  is  far  above  the  dynamic  range  capability  of  the  processor. 


V.  EVALUATION  OF  SYSTEM  PERFORMANCE 

In  the  previous  sections  we  developed  quantitative  estimates  of  the 
effects  of  various  component  errors  on  system  performance  --  namely  achievable 
cancellation  and  loop  stability.  In  this  Section  we  show  how  the  effects  of 
these  various  component  errors  reveal  themselves  in  measurements  ma^..  on  a 
real  system,  and  conversely,  how  one  can  tailor  the  system  measurements  to 
evaluate  the  errors  inherent  in  the  processor  or  antenna.  We  address  this 
subject  from  the  viewpoint  of  component  effects  on  cancellation.  Component 
effects  on  loop  stability  are  either  obvious  (e.g.,  unstable  loops)  or 
directly  measurable  by  observing  the  transient  behavior  of  the  loops  during 
adaption . 

An  experimental  adaptive  nulling  processor  was  constructed  and  evaluated 
in  order  to  discover  the  limits  on  performance  imposed  by  component  errors  in 
the  hardware.  A  detailed  description  of  the  processor  is  given  in  Ref.  9. 

Here  we  limit  our  discussion  to  the  measurement  techniques  employed  for 
evaluating  the  processor  cancellation  performance.  The  weight  circuits  were 
initially  the  PIN-diode  design  of  Fig.  9.  Following  that  a  set  of  weights 
based  on  transconductance  multipliers  was  installed  and  evaluated. 
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The  following  technique  was  used  to  evaluate  the  individual  weight 
circuits  before  they  were  installed  in  the  processor.  A  digital  feedback  loop 
was  implemented  using  an  HP  8505  network  analyzer  and  an  HP  9825  desktop 
computer.  A  CW  tone  at  band  center  is  applied  to  the  input  of  a  single  weight 
and  the  digital  feedback  loop  was  activated.  Closing  the  loop  causes  the 
weight  to  "turn  off."  At  this  point  the  weighted  path  exactly  cancels  the 
feedthrough  weighted  contribution  of  the  output  path  (see  Fig.  19).  The 
weight  I  and  Q  control  voltages  are  then  frozen,  and  the  input  is  swept  in 
frequency  across  the  band.  The  resultant  plot  of  output  power  vs  frequency 
for  two  of  the  weight  circuits  is  shown  in  Fig.  19.  It  is  clear  that  for 
wideband  systems,  the  feedthrough  can  contribute  significantly  to  the  total 
output  power  across  the  band.  These  weights  operated  at  an  IF  frequency  of 
120  MHz. 

Figure  19  illustrates  the  effect  of  weight  feedthrough  for  a  single 

channel.  In  a  processor  using  N  parallel  channels,  it  was  shown  in  Section 

2 

IV-C  that  the  cancellation  deteriorates  by  a  factor  of  N  (assuming  equal 

feedthrough  levels).  Measured  data  illustrating  this  are  presented  in 

Fig.  20.  Here,  one  channel  of  the  weight  sub-system  is  fixed  as  a  reference 

input,  and  the  number  of  cancelling  channels  is  varied  from  2  to  8  using  a 

switch  network  designed  into  the  test  system.  Broadband  noise  of  fixed 

bandwidth  is  applied  to  the  input  of  each  weight  through  an  N-way  divider. 

The  resulting  measured  output  spectra  for  various  numbers  of  cancellation 

2 

channels  is  illustrated  by  the  N  dependence.  For  the  same  measurement,  the 

effect  is  illustrated  in  a  different  manner  in  Fig.  21.  Here  the  measured 

increase  in  total  output  power  is  plotted  as  the  number  of  cancelling  channels 

2 

and  is  increased  compared  to  the  theoretical  N  behavior.  The  measurements 
indicate  that  three  of  the  weight  circuits  have  7.5  dB  more  feedthrough  than 
the  other  four. 

The  above  results  are  for  the  PIN-diode  weight  design  of  Fig.  9.  In  this 
circuit  the  feedthrough  is  dominated  by  the  isolation  and  time  delay  of  the 
quadrature  hybrid.  The  second  weight  design  we  evaluated  uses  the  principle 
of  transconductance  multipliers  and  is  illustrated  in  Fig.  22.  To  assess  the 


RELATIVE  OUTPUT  POWER  (0  dB  =  no  attenuation ) 


Fig.  19.  Weight  feedthrough  effects  on  cancellation  over  the  nulling 
band  using  a  simple  weight.  Results  for  two  different  channels  (4  and 
5)  are  plotted. 
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FIXED  REFERENCE  —  VARY  NUMBER 
OF  CANCELLING  CHANNEL 


CANCELLATION  vs  NUMBER  OF  CHANNELS 


Fig.  20.  Effects  of  weight  feedthrough  on  cancellation  as  the 
number  of  cancelling  channels  is  increased  from  two  to  eight. 


EASE  IN  OUTPUT  POWER(dB) 


Fig.  21.  Measured  increase  in  total  output  power  due  to  weight 
feedthrough  as  a  function  of  the  number  of  cancelling  channels. 
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multipliers.  Center  frequency  =  121.4  MHz. 


improvement  over  the  PIN-diode  weights,  we  use  the  eigenvalue  estimate  of  the 
achievable  cancellation  given  by  Eq.  (24).  With  fixed  I  and  Q  control  voltages 
the  amplitude  and  phase  response  of  each  channel  is  measured  separately  as  a 
function  of  frequency  with  a  network  analyzer.  The  correlation  matrix,  R,  is 
computed  from  the  data  parametrically  as  a  function  of  bandwidth,  BW.  The 
eigenvalues  of  R  are  then  computed  and  plotted  as  functions  of  BW.  The  result 
of  processing  the  data  in  this  way  is  illustrated  in  Fig.  23  for  both  the 
transconductance  multiplier  weights  and  for  the  PIN-diode  weights.  Recalling 
from  Eq.  (24)  that  C  “  s^/s^,  we  anticipate  a  15  dB  improvement  in  wideband 
cancellation  performance  with  the  transconductance  weights.  This  performance 
improvement  was  indeed  realized  when  the  weights  were  installed  in  the  processor. 

The  previous  discussion  has  indicated  several  techniques  useful  for  the 

evaluation  of  the  effects  of  weight  feedthrough.  We  have  shown  that  feed- 

2 

through  results  in  an  N  degradation  in  broadband  output  power  relative  to  a 
single  channel  even  if  the  weight  feedthrough  paths  are  identical  in  their 
frequency  response  (i.e.,  they  track  perfectly). 

We  now  consider  the  effects  of  channel  tracking  errors  which  arise  from 
the  fact  that  the  frequency  transfer  functions  of  each  channel  differ  slightly 
because  of  component  differences.  When  considering  the  effect  of  channel 
tracking  errors,  we  emphasize  that  it  is  the  total  channel  response  from  input 
to  output  which  must  track.  A  valuable  technique  for  separating  the  effects 
of  different  components  on  the  overall  tracking  performance  of  a  channel  is  to 
assemble  and  evaluate  the  processor  sequentially,  beginning  with  the  N-way 
combiner,  proceeding  to  the  weights,  then  the  IF  sections,  and  finally  the 
RF  front  ends,  until  the  total  processor  is  assembled.  At  each  stage  of 
assembly,  channel  tracking  measurements  with  data  reduction  similar  to  that 
discussed  above  are  made  to  determine  which  stages  dominate  the  cancellation 
performance . 
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4. 


EIGENVALUE  (dB) 


BANDWIDTH  (MHz) 


Fig.  23.  Eigenvalues  of  the  7-weight  channels  vs.  bandwidth  (MHz) 
for  both  the  pin-diode  weight  design  of  Fig.  9  and  the  transcon¬ 
ductance  multiplier  weight  design  of  Fig.  22.  Frequency  response 
of  each  channel  obtained  using  a  fixed  control  voltage  Vj  =  -1, 
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When  this  process  is  complete  it  will  be  found  in  general  that  one  or 
two  channels  are  measurably  worse  than  the  others.  Figure  24  illustrates 
the  effect  of  an  anomalously  bad  channel  on  the  eigenvalues  of  the  processed 
channel  tracking  data.  For  an  earth-coverage  quiescent  radiation  pattern 
and  a  planar  array  of  earth  coverage  elements,  only  a  single  channel  will  be 
"on"  in  the  unadapted  mode.  The  corresponding  beam  steering  vector  is 
V+  =  [0,  ...  0,1,0  ...]  where  the  1  appears  in  the  reference  (or  "on"l  channel. 
Choosing  a  dominant  error  channel  as  the  reference  degrades  the  performance 
significantly.  In  fact,  if  one  measures  the  cancellation  as  a  function  of 
reference  channel,  there  will  be  one  channel  which  yields  the  best  cancel¬ 
lation.  This  channel  should  then  be  designated  as  the  reference. 

Figure  25  illustrates  the  effect  for  the  transconductance  multiplier 
weight  circuits  of  the  processor.  It  shows  the  measured,  broadband  cancel¬ 
lation  as  a  function  of  reference  channel.  There  are  7  channels  in  all. 

Note  that  a  significant  improvement  in  cancellation  can  be  achieved  by 
selecting  the  reference  channel.  For  comparison.  Fig.  25  also  shows  the 
computed  cancellation  as  determined  from  the  eigenvalues  of  the  measured 
tracking  data  from  the  swept  frequency  measurement  discussed  previously. 

The  discrepancy  between  the  two  curves  indicates  that  tracking  errors  in 
these  particular  weights  are  a  function  of  weight  setting.  The  swept 
frequency  tracking  data  were  obtained  with  fixed  I  and  Q  control  drive 
voltages  to  each  weight  (in  their  minimum  attenuation  state) ,  whereas  the 
actual  cancellation  was  measured  with  the  weights  at  their  adapted,  steady- 
state  values.  Figure  25  applies  to  the  weight  circuits  alone. 

Figure  26  shows  the  equivalent  data  with  the  remaining  stages  of  each 
channel  added.  These  are  the  IF  amplifiers,  band-limiting  filters  and  finally 
the  RF  preamplifiers  and  mixers.  The  measured  and  computed  results  (again, 
computed  from  the  measured  tracking  data)  are  now  more  closely  in  agreement 
confirming  that  the  transconductance  multiplier  weights  are  not  the  dominant 
components  in  limiting  the  achievable  cancellation.  Of  the  remaining  IF  and 
RF  components,  there  was  in  fact  no  dominant  component.  Each  RF  and  IF  sub¬ 
section  contributed  in  roughly  equal  measure  to  the  cancellation  degradation. 
The  existence  of  an  "optimum"  reference  channel  is  also  evident  from  Fig.  26. 
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Fig.  24.  Effect  of  an  anomalously  bad  channel  on  the  cancellation 
performance  of  the  processor  when  the  bad  channel  is  used  as  a 
reference  channel. 
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CANCELLATION  (relalive  scale) 
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w-Z-1,-1  ON  REFERENCE  CHANNEL 
-10  OTHER  CHANNELS 


TRACKING 
VARIATION  WITH 
WEIGHT  SETTING 


REFERENCE  CHANNEL 

Fig.  25.  Measured  cancellation  vs.  channel  used  as  reference  for 
transconductance  multiplier  weights,  compared  to  the  computed  can¬ 
cellation  using  the  measured  frequency  response  for  each  channel 
with  control  voltage  =  (-1,  -1). 
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Fig.  26.  Measured  cancellation  vs.  channel  used  as  reference  with 
the  remaining  stages  (IF  through  UHF)  of  the  processor  added. 
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The  previous  results  have  illustrated  the  use  of  channel  tracking  meas¬ 
urements  for  determining  component  effects  on  achievable  cancellation.  Still 
another  presentation  exists  which  yields  the  effects  of  specific  component 
errors  on  performance,  and  this  technique  will  now  be  considered.  This  useful 
technique  can  be  explained  with  reference  to  the  output  power  vs  input  power 
characteristic  of  a  nulling  processor*  (Fig.  27).  Consider  first  an  ideal 
system  employing  a  hard  limiter  in  each  feedback  channel  operating  first  in 
the  quiescent  (without  nulling),  then  in  the  adapted  (with  nulling)  mode. 

Assume  that  a  test  signal  of  variable  bandwidth  and  power  is  applied  to  the 
processor  through  an  N-way  splitter.  This  simulates  broadside  incident  of  a 
single  source.  In  the  quiescent  mode,  as  to  a  given  channel  increases, 

P  increases  in  a  linear  relationship  to  P^.  In  the  adapted  mode,  a 
threshold  level,  P  ,  exists  where  the  nulling  processor  just  "senses"  the 
input  power  level.  At  this  threshold,  6  dB  of  cancellation  results.  As 
increases  further,  P  .  ideally  remains  fixed  at  a  constant  level  as  a  cor.se- 
quence  of  the  hard-limiting  action  in  the  feedback  loop.  The  achieved  cancel¬ 
lation  is  very  nearly  equal  to  the  difference  between  the  quiescent  and 
adapted  levels.  In  an  ideal  system,  therefore,  cancellation  increases  without 
limit  as  Pin  increases.  For  a  real  system,  however,  only  some  maximum  cancel¬ 
lation  is  achievable.  As  one  increases  P.  above  threshold  in  a  real  processor, 

m 

the  output  power  remains  constant  only  over  a  limited  range  of  P^  .  At  some 
value  of  ,  the  output  power  begins  to  increase  and  becomes  proportional  to 
Pin‘  repeating  this  measurement  as  each  section  of  the  processor  is  added, 
the  result  of  the  measurement  then  yields  directly  the  error  effect  that  is 
dominating  the  performance  degradation.  By  tailoring  the  input  appropriately, 
this  effect  can  usually  be  identified.  Representative  results  are  illustrated 
on  the  Figure,  not  necessarily  in  order  of  importance. 


*The  characteristic  curve  in  Fig.  27  is  applicable  to  processors  with  hard 
limiters  in  the  feedback  loops  (Ref.  7). 
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Figure  28  is  an  example  of  this  technique  as  applied  to  the  demonstration 
processor  using  the  PIN-diode  weight  circuits.  Several  cases  of  interest  are 
presented.  First  consider  the  earth  coverage  mode,  where  a  single  channel  is 
used  as  reference.  When  the  test  input  is  a  tone,  then  the  achievable  cancel¬ 
lation  is  limited  by  the  weight  noise  arising  from  the  baseband  circuitry 
(Section  IV-E-1).  When  a  wideband  test  signal  is  applied  to  the  processor 
simulating  broadside  incidence  on  the  array,  the  achievable  cancellation  is 
limited  either  by  channel  tracking  errors  or  weight  feedthrough.  Both  effects 
yield  the  same  trends  when  PQut  is  plotted  vs  P^  .  Therefore  this  technique 
alone  cannot  distinguish  between  the  two.  The  fact  that  it  was  weight  feed¬ 
through  from  the  PIN-diode  weights  that  dominated  these  results  was  determined 
2 

by  the  N  dependence  discussed  previously.  Another  interesting  fact  is  that 
when  the  steering  vector  changes  from  earth  coverage  to  a  maximum  gain  beam 
mode,  the  quiescent  output  power  increases  accordingly,  because  the  array  gain 
increases.  However  (as  discussed  in  Section  I I 1-3) ,  the  adapted  output  in  the 
feedthrough  dominated  region  is  independent  of  the  steering  vector.  This 
results  in  greater  cancellation  for  one  steering  vector  than  for  another. 

This  is  an  example  of  the  ambiguity  inherent  in  using  cancellation  level  as  a 
means  of  specifying  system  performance. 

Figure  29  illustrates  similar  results  for  a  case  where  the  test  signal 
was  time  delayed  to  each  channel  in  a  way  to  simulate  an  8°  angle-of-arrival 
on  a  specified  array  geometry  to  emphasize  the  effect  of  array  dispersion. 
Comparison  of  the  measurement  with  computer  simulation  of  the  array  bandwidth 
characteristics  yields  results  which  confirm  that  array  dispersion  is  the 
dominant  factor  in  this  case. 

The  previous  discussion  has  centered  on  evaluating  the  effects  of  compo¬ 
nent  errors  (i.e.,  deviations  from  ideal  behavior)  on  performance.  We  now 
consider  the  effect  of  hard-limiting  in  the  feedback  loop.  While  not  an  error 
source  strictly  speaking,  hard  limiters  can  have  a  harmful  effect  on  perform¬ 
ance  and  for  that  reason  are  mentioned  briefly  here.  For  a  single  interfer¬ 
ence  source,  it  is  relatively  simple  to  show  that  hard- limiting  reduces  the 
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OUTPUT  POWER  (relative  scale) 


Fig.  28.  Measured  cancellation  trends  —  output  power  vs.  input  power 
for  selected  inputs  and  beam- steering  parameters. 


OUTPUT  POWER  -  RELATIVE  SCALE 


/ 


INPUT  POWER -RELATIVE  SCALE (dB) 


Fig.  29.  Measured  cancellation  trends  -  output  power  vs.  input  power 
for  input  simulating  an  8°  angle  of  incidence  using  the  7-element 
double  triangle  configuration  in  the  insert  of  Fig.  17. 
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circuit  dynamic  range  requirements  in  the  feedback  loops.  For  example,  if  the 

processor  is  designed  to  operate  over  a  40  dB  dynamic  range  of  input  power, 

the  correlator  mixer  and  subsequent  circuits  must  operate  over  an  80  dB  range 

of  power  if  no  hard  limiter  is  employed.  With  hard  limiting,  the  dynamic 

range  requirement  on  the  correlator  is  reduced  to  40  dB.  However,  if  there 

are  two  interference  sources  which  are  significantly  spread  in  power  level, 

they  give  rise  to  two  widely  spaced  eigenvalues  in  the  signal  correlation 

C71 

matrix.  It  can  be  shown  that  the  effect  of  the  hard-limiter  is  to 
suppress  the  smaller  eigenvalue  (see  Table  1),  so  that  it  is  essentially  not 
sensed  by  the  processor.  The  degree  to  which  this  happens  depends  on  the 
degree  to  which  the  larger  eigenvalue  is  above  threshold.  For  this  measure¬ 
ment,  the  eigenvalue  s1  (controlled  by  an  incident  source  at  a  given  angle 
of  arrival  having  power  pp  is  fixed  relative  to  threshold.  The  second  eigen¬ 
value  s^  (controlled  by  a  second  source  eminating  from  a  different  angle  of 
arrival  having  power  P2)  is  varied,  and  the  output  power  is  measured  as  s2  is 
varied.  The  results  are  plotted  in  Fig.  30  for  a  system  having  a  specific 
dynamic  range.  The  region  over  which  the  small  eigenvalue  is  not  sensed  is 
determined  by  Sj/Pj.,  the  ratio  of  the  larger  eigenvalue  relative  to  threshold. 
By  examining  the  measured  results,  it  becomes  clear  that  the  system  cancella¬ 
tion  taking  into  account  both  sources  is  limited  to  •v  ✓s1/Pt  .  Hence,  if 
Sj/Pt  is  the  desired  cancellation  level,  then  the  dynamic  range  required  of 
the  baseband  circuitry  is  s1/Pt  whether  or  not  a  hard-limiter  is  used.  A 
small  eigenvalue  might  arise  from  a  second  incident  source  of  low  power  level 
(as  simulated  above) ,  or  it  may  result  from  a  single  wideband  source  incident 
on  a  dispersive  antenna.  In  this  second  case  the  hard- limiter  essentially 
prohibits  the  processor  from  compensating  for  array  dispersion  by  using  up 
more  than  one  degree  of  freedom  to  form  a  broadband  null.  In  summary,  a 
hard-limiter  does  not  reduce  the  circuit  dynamic  range  requirements  when  two 
sources  of  widely  differing  power  levels  are  present,  and  requires  a  broad¬ 
band  antenna  geometry  when  deep  cancellation  of  a  single  source  is  desired. 

We  conclude  that  processors  which  must  null  multiple  sources  over  a  wide 
dynamic  range  of  power  should  not  use  hard  limiters  in  the  feedback  loops. 
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OUTPUT  POWER -RELATIVE  SCALE 


W  «i\»  T* 


Fig.  30.  Measured  hard-limiter  suppression  of  the  second  eigenvalue 
generated  by  a  weaker  interference  source. 
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APPENDIX 


In  this  Appendix  we  estimate  the  eigenvalue  ratio  S2/s^  for  a  phased 
array  antenna  employing  N  elements  each  having  element  gain  Gg.  Let  a  wave- 
front  be  incident  on  the  array  from  angles  0,  4>.  The  wavefront  is  normalized 
to  unity  amplitude.  The  correlation  matrix  as  defined  at  the  antenna  output 
ports  is  given  by  (for  broadband  noise  incidence) 

-T  ) 

R  =  <Ge  e  q  >  (Al) 

where  denotes  the  frequency  average 

u>o+BW/2 

<*>  -  1W  /  dw  (A2) 

a)  -BW/2 
o 

th  * 

and  is  the  time  delay  for  the  wavefront  to  reach  the  k—  element  relative 
to  an  arbitrary  reference  position  on  the  array: 

Tk  =  ~  sinejxk  cos$  +  y^  sin^J  (A3) 

D  2  2 

and  the  x^,  y^  are  the  element  position  normalized  to  (x^  +  yk  ±  1).  The 
eigenvalue  Sj  is  simply  given  by  NGg,  since  Sj  corresponds  to  applying  a 
set  of  weights  which  points  a  maximum  beam  toward  the  interferfence  source 
having  maximum  array  gain  NGg.  In  order  to  estimate  S2,  we  expand  R  in  the 
form 

R  =  Rq  +  AR  (A4) 


where  Rq  is  the  correlation  matrix  assuming  zero  bandwidth,  and  AR  represents 
the  perturbation  to  R^  due  to  non- zero  bandwidth  effects.  Clearly,  then, 
from  Eq.  (Al), 


Gee 


-HCVV 


(A5) 


Al 
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1 


To  obtain  AR,  we  expand  R  in  the  form 


-jVW 

R.  =  Ga<e  °  k  q  e  k  q  >  a 

=k,q  eN  ' 

‘ja>o(W  ^“o(W 

G  e  q  -  j  G  <Aw>  (T. -T  )  e  °  q 


‘k  V 


-Ge  <r>(VTq)^e 


2  ^WV 


The  term  involving  Aw  averages  to  zero,  so  that  from  (A6)  we  obtain 

(w  FBW)  2  -iWV 

AR  =  -  G  — ~ —  (T. -T  )  e  q 

=  e  24  v  k  q 

The  eigenvalue  s2  can  then  be  estimated  according  to: 

1* 

Max  £  •  R  •  £ 

S2  =  i  - t“ - 

e  •  e 


(A7) 


(A8) 


where  only  vectors  £  such  that  £  •  £.  =  0  are  used  in  the  maximization. 
Since  £  •  R  «  0,  (£^  is  the  first  eigenvector  of  Rq  and  R) ,  then 

Max  (w  FBW)2  -jw  T,  )  +ju>  T  ■> 

*2  ■  •  Ge  -V-  £S>k*  •  °  “  <•,  •  °  <Vy2 

k  q  M  M 

We  note  from  (A9)  the  special  case  (broadside  incidence)  where 

s2  =  0.  Furthermore,  since 

■jw  T. 

L.  J  o  k 
A 


(A9) 


6k  6 


=  0 


(A10) 


t  2 

(i.e.,  £  •  Rq  •  e  =0),  the  factor  (T^-T  )  in  (A9)  essentially  looks  like 
a  correlated  tracking  error  between  channels  which  increases  as  (D/X  •  sin©)*-. 
Furthermore,  Eq.  (A10)  is  independent  of  D/X,  as  ef  simply  negates  the  D/X 
-juoTk 

phase  term  e  K  so  that  the  sum  in  (A10)  adds  to  zero.  Thus  the  only 
term  inside  the  sum  in  (A9)  which  depends  on  D/X  is  the  factor 
Noting  Eq.  (A3),  we  can  factor  D/X  sin©  from  inside  the  sum  to  rewrite  (A9) 
in  the  form: 
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